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1.  Purpose.  The  purpose  of  this  memorandum  is  to  provide  documentation  of  research 
for  the  Army  G3/5/7  Munitions  Management  Division  and  Center  for  Army  Analy¬ 
sis  (CAA)  by  the  TRADOC  Analysis  Center,  Monterey  (TRAC-MTRY).  The  focus 
of  the  research  is  to  improve  the  current  methodology  employed  by  CAA  in  model¬ 
ing  precision-guided  munitions  effects  into  the  Joint  Integrated  Contingency  Model 
(JICM)  scenarios  to  enhance  munitions  requirements  estimate  generation,  especially 
for  high-cost  low-density  munitions  like  Excalibur. 

2.  Background.  The  Army  Munitions  Management  Division  develops  ammunition  re¬ 
quirements  in  support  of  Army  programming  and  the  Army  budget.  When,  as  usually 
happens,  there  are  not  sufficient  resources  to  fill  100%  of  requirements,  this  office  also 
offers  analysis  on  the  risks  associated  with  not  filling  the  full  requirement  by  munition. 
Currently,  CAA  supports  the  Munitions  Management  office  by  using  JICM  to  do  Army 
level  analysis  to  support  requirements  and  risks  associated  with  tradeoffs  in  filling  re¬ 
quirements.  CAA  uses  effects  data  provided  by  the  Army  Material  Systems  Analysis 
Activity  (AMSAA)  and  tactics,  techniques,  and  procedures  (TTPs)  and  related  data 
from  subject  matter  experts  at  the  Fires  Center  of  Excellence  (FCoE)  to  inform  the 
development  of  requirements  for  precision  guided  munitions  (PGMs).  Given  the  criti¬ 
cal  role  that  PGMs  play  in  modern  Joint  warfare,  it  is  of  great  importance  to  ensure 
that  CAA’s  modeling  represents  PGM  usage  and  effectiveness  in  an  operationally  ap¬ 
propriate  manner  to  inform  risk  analysis  for  not  hlling  precision  munitions  to  100%  of 
requirements. 

3.  Methodology.  The  research  team  used  two  approaches  to  solve  the  JICM  precision- 
munitions  effects  modeling  problem:  an  analytic  model  for  the  Carleton  damage  func¬ 
tion,  and  a  response  surface  model  to  predict  the  probability  of  damage,  wounded 
personnel,  suppression  effects.  Appendix  [A]  shows  constraints,  limitations,  and  as¬ 
sumptions  used  for  our  research  methodology. 

(a)  The  first  approach  refined  the  analytical  solution  to  the  Carleton  damage  func¬ 
tion  to  compute  the  damage  probability  for  multiple  munitions  fired  at  a  unitary 
target  and  area  targets.  Wang  et  al.  in  detail  the  derivation  of  their  equations 
(See  Appendix  0.  First,  Wang  et  al.  determined  the  requisite  parameters  to 
determine  an  exact  damage  probability  using  the  Carleton  damage  function  of  M 
weapons  fired  at  a  unitary  target.  Next,  they  utilized  this  solution  to  optimize 
the  aim-point  distribution  that  provides  the  greatest  probability  of  damage  to  the 
target 0  The  results  of  this  research  will  allow  the  JICM  developers  to  improve 
their  modeling  methodology  which  does  not  evaluate  the  optimal  probability  of 
multiple  weapons  systems  fired  at  a  unitary  target.  After  determining  the  optimal 
damage  probability  against  a  unitary  target,  Wang  et  al.  identified  the  optimal 
probability  against  an  area  target]^]  Essentially,  this  research  extends  their  previ- 

1Hongyun  Wang  et  al.  “Explicit  Exact  Solution  of  Damage  Probability  for  Multiple  Weapons  against  a 
Unitary  Target”.  In:  American  Journal  of  Operations  Research  6.06  (2016a),  p.  450. 

2Hongyun  Wang  et  al.  Average  damage  caused  by  multiple  weapons  against  an  area  target  of  normally 
distributed  elements.  Tech.  rep.  TRAC-Monterey,  2016b. 
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ous  work  to  multiple  target  engagements  and  provides  added  capability  for  the 
JICM  developers  to  implement  in  their  combat  model. 

(b)  Although,  the  analytic  methodology  provided  added  benefits  to  the  JICM  model 
by  improving  the  number  of  munitions  expended  in  the  model,  it  still  does  not 
ensure  the  fidelity  of  Battlefield  effects  required  by  the  Army  G-3/5/7  munitions 
management  personnel.  The  response  surface  modeling  methodology  addresses 
this  issue  by  computing  a  predicted  probability  of  damage,  wounded  personnel, 
and  suppression  based  on  the  current  inputs  used  in  the  JICM  model.  (See  Ap¬ 
pendix  [Dj.  To  start,  Ahner  and  McCarthy  computed  a  damage  probability  using 
the  Klopcic  damage  function,  which  is  a  combination  of  the  Carleton  damage 
function  and  the  cookie  cutter  function^]  Next,  they  used  this  damage  function 
and  the  JICM  input  parameters  to  build  a  design  of  experiments  and  a  meta¬ 
model  to  compute  a  predicted  probability  of  damage  based  on  the  munition  type 
and  target  distance.  They  also  developed  a  similar  methodology  to  compute  a 
predicted  probability  of  suppression  and  wounded  personnel j^]  Finally,  Ahner  and 
McCarthy  developed  a  risk  value  to  model  a  commander’s  risk  tolerance  for  col¬ 
lateral  damage.  Using  this  parameter  will  give  a  decision  maker  insight  into  the 
amount  of  certainty  they  can  have  of  damaging,  wounding,  or  suppressing  a  partic¬ 
ular  target  which  is  currently  not  a  capability  within  the  JICM  model.  Another 
benefit  of  their  methodology  is  that  the  CAA  JCIM  developers  can  implement 
this  algorithm  into  their  current  model  or  use  a  separately  developed  spreadsheet 
model  as  a  supplemental  analysis  tool. 

4.  Results.  The  project  team  briefed  the  final  results  to  the  sponsor  on  27  October  2016 
and  the  slides  and  notes  from  this  briefing  are  contained  in  Appendix  |B| 


3Darryl  Ahner  and  Andrew  McCarthy.  Algorithm  Development  for  the  Combat  Sample  Generator 
(COSAGE)  Model.  Tech.  rep.  Air  Force  Institute  of  Technology,  2016. 
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Appendix  A 
Study  Plan 


Problem  Statement 

To  improve  the  current  methodology  employed  by  the  Center  for  Army  Analysis  (CAA) 
in  modeling  precision-guided  munitions  effects  into  JICM  scenarios  in  order  to  enhance 
munitions  requirements  estimate  generation,  especially  for  high-cost  low-density  munitions 
like  Excalibur. 


Project  Team 

Sponsor  Agency:  Robert  Grubbs 

Munitions  Management  Division,  HQDA  G-3/5/7 
robert .  a.  grubbs.  civ@mail .  mil 


Project  Stakeholder:  David  Knudson 

Campaign  Analysis  Division  Chief 

Center  for  Army  Analysis,  Fort  Belvoir,  VA  22060 

david.b.knudson.civ@mail.mil 


TRAC  Lead: 


NPS  Faculty: 


AFIT  Faculty: 


Cardy  Moten  III 
MAJ,  LG/FA49 

TRADOC  Analysis  Center  -  Monterey 
cardy.moten3.mil@mail.mil 

Dr.  Hong  Zhou 

Applied  Mathematics  Department 

Naval  Postgraduate  School,  Monterey,  CA 

hzhou@nps.edu 

Dr.  Darryl  Ahner 

Operations  Research  Department 

Air  Force  Institute  of  Technology,  Dayton,  OH 

darryl .  ahner  @  aht .  edu 
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Constraints,  Limitations,  &  Assumptions 

•  Constraint 

—  The  project  completion  date  is  no  later  than  30  September  2016. 

•  Limitation 

—  The  project  will  focus  only  on  precision-guided  munitions  estimates  and  effects. 

•  Assumption 

—  THe  project  methodology  can  be  implemented  for  other  types  of  munitions  within 

JICM. 


Methodology 
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Timeline 

06  JAN  16 

15  JAN  16 
18  JAN  16 
9  MAR  16 
17  MAY  16 
27  SEP  16 
27  OCT  16 

16  DEC  16 


IPR  #1:  Restated  problem  statement  approval. 
Complete  literature  review. 

Develop  modeling  methodology. 

IPR  #2:  Brief  initial  findings  to  CAA. 

IPR  #3:  Brief  initial  findings  to  G3  Munitions. 
Final  IPR  to  CAA. 

Final  IPR  to  G3  munitions. 

Deliver  project  documentation. 
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Appendix  B 
Progress  Report 


The  final  IPR,  presented  to  the  sponsor  on  27  October  2016,  for  this  phase  of  the  project  is 
on  the  following  pages. 
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In  the  current  process,  JICM  produces  an  estimate  of  Excalibur  rounds  fired.  This  estimate  is 
manually  mapped  to  FIRESIM  results.  The  problem  with  this  method  is  that  FIRESIM  and 
JICM  do  not  adjudicate  their  effects  at  the  same  resolution. 

Since  FIRESIM  is  a  higher  resolution  model  than  COSAGE,  the  current  process  requires  more 
subjective  analysis  of  the  results. 
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[E3I  CAA  needs  this  capability  to  provide  robust  munitions 
S3  estimates  to  the  Army  G3  Munitions  Division. _ 


Munitions 

Estimate 

r 

t 

Updated 

Constraints 

UNCLASSIFIED 


This  slide  makes  a  simple  depiction  of  why  this  problem  is  important.  In  essence,  CAA  sends 
an  estimation  of  munitions  requirements  to  the  Army  G3  Munitions  office  for  their  follow-on 
analysis. 


The  Army  G3  munitions  office  provides  feedback  to  CAA  in  the  form  of  updated  munitions 
constraints  for  CAA  to  model  and  then  provide  updated  munitions  estimates. 


Since  the  estimates  for  precision-guided  munitions  are  not  adjudicated  in  enough  detail  to 
provide  the  Army  G3  munitions  division  an  assessment  of  the  relative  risk  to  mission,  the 
Army  G3  munitions  office  is  having  difficulty  in  providing  adequate  constraints  model  for 
CAA  to  improve  the  munitions  estimates. 


Thus,  our  proposed  model  should  provide  a  potential  solution  to  this  problem. 
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[no] 


COSAGE  needs  an  improved  effects  modeling 
methodology  for  precision-guided  munitions. 
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The  current  COSAGE  methodology  is  as  follows. 

-  Select  1  piece  of  critical  equipment  to  target  at  random 

-  Based  on  range  from  target  a  lethal  area  and  associated  round  errors  are  determined  for  the 
piece  of  equipment  selected  at  random 

-  Calculation  of  PK  for  1  round 

Center  point  of  piece  of  equipment  is  the  aim  point  (0,0) 

Target  size  is  based  on  the  largest  dimension  of  the  target  equipment  size 
Use  a  D0  value  of  1 .0 

Use  the  Carlton  Damage  Function  from  ME  to  determine  a  PK 

-  Draw  a  random  number  to  determine  if  the  piece  of  equipment  was  killed  or  not 

There  is  no  Collateral  targets  under  consideration 

To  improve  this  methodology,  we  are  going  to  explore  more  robust  damage  functions,  and 
improve  the  target  location  algorithms. 
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[Ell  ^se  our  m°dels  to  guide  your  development  of  a 
"  precision-guided  munitions  effects  model  in  COASE. 


- Depedent  error 

—  —  —  Inoepenflent  error 


Analytical  Effects  Model  Response  Surface  Effects  Model 
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The  bottom  line  is  that  we  feel  that  using  our  proposed  model  will  help  CAA  to  fill  its 
current  gap  in  adjudicating  precision-munitions  effects  in  COSAGE. 
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[Ell  At  endstate,  CAA  will  have  the  methods  to  develop  an  „ 

msm  improved  effects  model  for  precision-guided  munitions. 


JICM 

Model 


Produce 


Munitions  Estimate 


Battlefield  Effects: 

•  Damaged  Equipment 

•  Killed  personnel 

•  Wounded  personnel 
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The  intent  of  our  modeling  methodology  is  to  not  only  provide  information  on  how  many 
rounds  were  fired,  but  also  other  battlefield  effects  such  as:  damaged  equipment,  personnel 
killed,  and  wounded  personnel. 
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[E31  ^ur  m°deling  approach  is  broad  in  perspective  to 
[02=0]  ensure  we  analyze  the  most  important  issues. 

Conduct  Literature  Review 

Review  JICM  literature 

Review  precision-guided  munitions  literature 

Revised  problem  statement 


Develop  munitions  effects  model 

Develop  analytical  model. 

Develop  response  surface  model. 

What  is  the  recommended  method  to  represent  the  effects  of  the  impact  of  an 
Excalibur  munition  in  a  combat  model  that  has  a  resolution  above  aggregate  level? 
What  are  the  effects  for  a  point  target? 

What  are  the  effects  for  an  area  target? 

How  can  we  represent  effects  for  damage,  suppression,  and  wounded  casualties? 
Psuedo  code  and  algorithm  for  model  development. 

Technical  memorandum. 
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[C3I  We  use(^  the  following  constraint,  limitation,  and  an 
33  assumption  to  scope  our  problem. _ 

Constraint. 

The  project  completion  date  is  no  later  than  30  September  2016. 

Limitation. 

The  project  will  focus  only  on  precision-guided  munitions  estimates  and  effects. 

Assumption. 

The  project  methodology  can  be  implemented  for  other  types  of  munitions  within 
JICM. 


October  2016 _ UNCLASSIFIED _ Preci»»on~guKfcd  Munrton»tWect» 
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o 

Analytical  Model  Results 


- Depedent  ecror  I 

-  -  -Inoepenoerrt  error 
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Professor  of  Applied 
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ts  We  started  our  aPProach  by  analyzing  multiple 
BBS  weapons  against  a  unitary  target. _ 

M  weapons  =  either  multiple  weapons  or  multiple  shots  from  one  weapon 

Effect  of  each  weapon  is  modeled  by  Carleton  damage  function 

Dependent  error  =  one  error  deviating  impact  points  of  all  M  weapons 

Independent  errors  =  M  independent  errors,  one  affecting  each  of  M  weapons 


- Depedent  error 

- Independent  error 


Impact  point  of  weapon  j 


Wj 

Impact 

point 


r,  +  Y  + 


Aiming  Dependent 
point  error 


Independent 

error 
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For  this  case  we  consider  the  problem  of  shooting  multiple  weapons  against  a  single  target. 
We  model  the  effects  of  the  munitions  using  the  Carleton  damage  function.  We  also  will  model 
the  effects  by  analyzing  both  dependent  error  (aiming  error)  and  independent  error  (ballistics 
errors). 

In  short,  the  impact  point  for  the  round  is  a  combination  of  the  aiming  point  and  the  effects 
from  both  the  dependent  and  independent  error. 
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[E3I  Next’  we  computed  the  kill  probability  in  the  presence 
SB  of  dependent  and  independent  errors. _ 

Many  realizations  of  firing  M  weapons: 


Kill  probability  =  average  over  all  realizations 

Exact  solution  for  the  kill  probability: 

fast  and  accurate  calculation  of  kill  probability 

making  it  possible  to  optimize  aiming  points  of  M  weapons 

UNCLASSIFIED 


We  derived  our  kill  probability  solution  from  the  Carleton  damage  function  and  will  use 
this  solution  to  examine  the  kill  probability  corresponding  to  various  distributions  of  the 
aimpoints  of  M  weapons. 
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E3I  Essential|y>  spreading  out  aiming  points  is  effective  in 
88  countering  dependent  error. _ 


Survival  probability  =  1  -  (kill  probability) 


Survival  probability  does  not  decay 
exponentially  with  number  of  weapons. 
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The  left  panel  shows  a  comparison  of  the  case  of  all  M  weapons  aimed  at  (0,0)  vs  the  case 
of  using  the  optimal  distribution  of  aim  positions  in  the  decay  of  l-pkill.  With  the  optimal 
distribution  of  aim  points  for  M  weapons,  we  expect  that  l-pkill(M)  decays  faster  than  in  the 
case  of  aiming  all  M  weapons  at  the  same  position,  which  the  left  side  chart  confirms. 


In  the  right  panel,  we  plot  log(l-pkill)  vs.  M.  It  is  clear  that  in  the  presence  of  dependent 
error,  the  survival  probability  decreases  slower  than  the  geometric  decay. 
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[E3J  Next,  we  ana^yzed  an  area  target  with  randomly 
SH  distributed  elements. _ 

Many  realizations  of  target  elements  distribution 


s:  radius  of  area  target 

=  standard  deviation  of  target  elements 
distribution 

Average  damage  fraction: 

^damage  (area  target,  M  weapons)  =  Average  over  aN  realizations 

_ UNCLASSIFIED _ 


This  chart  shows  an  example  of  the  type  of  problem  we  analyzed.  In  each  graphic,  there  are 
12  targets  normally  distributed  around  the  center  of  the  target  located  at  (0,0).  Similar  to  the 
point  target  problem,  the  solution  for  the  damage  to  these  area  targets  are  an  average  over  all 
realizations  of  the  targets. 
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Next,  we  analyzed  a  constrained  optimization  over 
SssS  various  patterns  of  aiming  points. _ 


Patterns  of  B-family:  uniform  in  polar  (geometric)  angle  along  an  ellipse 
B1  B2  B3 


3  parameters  3  parameters  4  parameters 

In  contrast,  unconstrained  optimization  of  12  weapons  has  24  parameters! 


Patterns  of  A-family: 

uniform  in  parametric  angle  along  an  ellipse 
A1  A2  A3 


(j- 1) 

9=6+— - -2n , 

1  M 


\x  =r  cos 9 

\  i  *  i 

\y .  =  r  sin# 

[■'  1  y  i 
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The  goal  of  our  method  is  to  find  simple  and  efficient  “empirical”  methods  for  calculating 
nearly  optimal  aiming  positions.  This  approach  greatly  simplifies  the  numerical  complexity  of 
finding  the  optimal  aiming  points  at  the  price  of  obtaining  an  approximate  optimum.  We 
considered  6  constrained  patterns.  The  patterns  of  B-family  aimpoints  are  M  points  on  an 
ellipse  that  are  uniform  in  polar  angle.  The  pattern  of  A-family  aimpoints  analyzed  M  points 
along  an  ellipse  that  are  uniform  in  parametric  angle. 


Pattern  1  for  both  groups,  depict  M  points  uniform  in  either  parameter  or  polar  angle. 
Pattern  two  depicts  one  target  located  at  (0,0)  and  the  remaining  (M-l)  targets  uniform  in 
either  parameter  or  polar  angle.  Pattern  3  depicts  two  targets  along  the  x-axis  and  the 
remaining  (M-2)  targets  uniform  in  either  parameter  or  polar  angle. 


Patterns  1  and  2  have  three  parameters:  the  parameter/polar  weapon  angle,  effective  radius 
of  the  ellipse,  and  the  aspect  ration  of  the  ellipse.  Patter  3  also  contains  these  three  parameters 
plus  a  fourth  parameter  for  the  aiming  point  for  the  x-axis  targets. 
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The  left  panel  plots  the  optimal  target  damage  fraction  as  a  function  of  several  values  for 
the  target  radius  (s)  and  a  fixed  value  for  the  number  of  weapons  M. 


A  practical  question  regarding  resource  allocation  is  the  following.  Given  the  radius  of  area 
target  (s),  what  is  the  minimum  number  of  weapons  needed  to  achieve  a  given  threshold  of 
damage  fraction?  This  question  is  answered  in  the  chart  on  the  right.  This  graphic  shows  that 
for  any  given  threshold  of  damage  fraction,  the  minimum  number  of  weapons  needed  is  an 
increasing  function  of  the  area  target  radius  (i.e.,  larger  area  target  requires  larger  number  of 
weapons),  which  again  is  reasonable  and  consistent  with  our  intuition. 
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Thus,  aiming  points  predicted  by  scaling  are  very 
efficient  and  accurate  for  area  targets. 


We  need  to  calculate  only  one  set  of  optimal  aiming  points 


{r(s),j  =  1,2 . m}  =^.{r.(s0),  7  =  1,2, . m} 

aiming  points  predicted 
by  scaling  law 
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For  an  area  target  of  radius  s,  we  simply  calculate/predict  a  set  of  nearly  optimal  aiming 
points  using  the  scaling  law  shown  above.  We  evaluate  the  performance  of  this  efficient 
method  by  examining  the  damage  fraction  values  achieved  by  these  sets  of  nearly  optimal 
aiming  points.  Specifically,  for  each  area  target,  we  calculate  the  damage  fraction  values 
corresponding  to  three  sets  of  aiming  points: 

•  Aiming  points  calculated  in  the  unconstrained  optimization. 

•  Aiming  points  calculated  using  the  scaling  law. 

•  All  aiming  points  =  (0,0). 


The  charts  in  the  slide  compare  the  damage  fraction  values  caused  by  the  3  sets  of  aiming 
points  described  above  for  M  =  6  weapons  (left  panel)  and  for  M  =  10  weapons  (right  panel). 
The  damage  fraction  achieved  by  the  set  of  nearly  optimal  aiming  points  calculated  using 
scaling  is  indistinguishable  from  that  achieved  in  the  unconstrained  optimization  (true 
optimum)  while  the  damage  fraction  corresponding  to  all  weapons  aiming  at  (0;  0)  is  much 
lower.  Therefore,  we  conclude  that  the  scaling  law  is  an  efficient  and  accurate  method  for 
calculating  a  set  of  nearly  optimal  aiming  points. 
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[E3]  Goals  are  to  find  an  approximate  solution  to  the  kill  probability  for  A  ». 
[ami]  a  known  set  of  conditions  and  expected  munition  requirements 


?.. 


Cookie  Cutter  Function  is  simple  to  implement,  but 
underestimates  effects  of  munitions  on  distant  targets 


Mtoctton  |n  |  .«  Rar.yc  |n, 

Carleton  Damage  Function  tends  to  underestimate 
damage  in  close  proximity  to  the  detonation  point. 


Does  an  excellent  job  approximating  damage  far  from 
detonation  point 


Klopcic  Hybrid  Function  uses  best  of  both  functions 
Source:  Klopcic,  1990 
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Carlton  Damage  Functions  calculates  PK  as  a  function  of 
distance  from  blast 
munition  impact  angle 
direct  hit  effectiveness 
&  weapon  lethal  area 

But  the  Carlton  damage  function  tends  to  underestimate  damage  in  close  proximity  to  the 
detonation  point  but  does  an  excellent  job  approximating  damage  far  from  detonation  point 

In  contrast,  the  cookie  cutter  function  estimates  damage  a  constant  damage  for  close 
proximity  points  and  PK  =  0  for  points  farther  away. 

This  becomes  especially  important  when  analyzing  fratricide  and  collateral  damage, 
damage  far  away  from  detonation  is  an  important  consideration,  which  is  more  appropriate 
modeled  by  the  Carlton  Damage  Function. 

We  therefore  combine  them  taking  the  best  characteristics  from  each. 
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The  damage  function  accounts  for  typical  inputs  to  produce  the  effects  of  one  munition.  For 
multiple  munitions,  the  effects  are  considered  as  independent  events  and  is  applied  repetitively. 
Terrain  affects  are  accounted  for  through  the  Lethal  Area  and  Core  Lethal  Area  inputs. 


Target  location  area  is  the  only  uncertain  factor  entering  into  our  equation.  We  treat  it  as  a 
random  variable  that  is  distributed  uniformly  from  0  to  its  maximum  value  in  all  directions. 
This  requires  a  simulation  to  determine  the  average  effect  since  our  damage  function  is 
nonlinear. 
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Munitions  affect  targets  and  nontargeted  objects  in 
various  ways 
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Precision  munitions  can  affect  any  object  in  their  area.  A  target  or  other  object,  such  as  a 
civilian,  can  be  affected  over  the  entire  spectrum  from  unaffected  to  killed.  An  enemy  who  is 
suppressed  is  unable  to  fire  and  there  may  be  a  time  dependent  component  to  this  effect. 

An  enemy  who  is  wounded  may  perform  at  less  than  peak  effectiveness,  but  may  or  may 
not  be  suppressed 

An  enemy  who  is  killed  is  eliminated  from  the  scenario 

The  probability  of  being  wounded  or  suppressed  is  modeled  similarly  to  P(K),  using  a 
scaler  to  determine  the  size  of  the  “wounding  lethal  area” 
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gg  Inputs  are  varied  to  ensure  coverage  of  all  operational 
xsrm  possibilities  when  constructing  the  simulation  runs 


Input  Variable 

Minimum 

Maximum 

Deflection  Distance 

0  meters 

50  meters 

Range  Distance 

0  m 

50  m 

D0 

0.05 

1.0 

Lethal  Area  (AL) 

100  m2 

1500  m2 

Core  Lethal  Area 

10%  of  Al 

60%  of  Al 

Terrain  factor 

0.5 

1.0 

ImpactAngle 

O 

O 

CO 

CD 

O 

o 

Target  Location  Error 

2m 

91  m 

Wound  LA  factor 

1.1 

2.5 

Suppression  LA  factor 

1.1 

2.5 
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ggj  Nearly  orthogonal  Latin  hypercube  is  used  on  inputs  to 
SH  create  response  surface  over  entire  space  of  inputs. 


Source:  Cioppa  and  Lucas,  2007 
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The  experimental  design  is  a  Nearly  Orthogonal  Latin  Hypercube  (NOLH),  which  is  a 
space  covering  design  that  ensures  the  correlation  between  input  variables  will  be  negligible 
which  has  better  space-filling  properties  than  completely  randomized  variable  selection 

Variable  correlation  is  restricted  to  being  in  the  range  (-.03,  .03) 

A  minimum  and  maximum  value  must  be  set  for  each  input  variable 
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Notice  that  we  end  up  modeling  the  square  root  of  Pk  in  order  to  satisfy  fitting  a  linear 
model.  RA2,  the  measure  of  model  fit  want  close  to  one,  changes  from  0.66  to  0.79  with  normal 
distribution  of  errors  around  the  solid  red  line. 

Unfortunately,  there  are  still  too  many  small  values  of  PK  that  affect  fit  resulting  in  poor 
predicted  larger  PK  values  that  we  care  about  most.  We  therefore  decided  to  use  regression 
trees  to  construct  the  model  piece  wise. 

Models  were  developed  for  different  maximum  values  of  TLE,  Lethal  Area,  and  distance 
from  the  origin  (d  and  r) 

Max  TLEs  of  15m,  50m,  and  91m 

Max  Lethal  Areas  of  500m2,  1000m2,  and  1500m2 

Max  d  and  r  of  10m,  15m,  20m,  25m,  and  50m 
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The  regression  tree  analysis  showed  that  Deflection  and  Range  were  the  two  factors  that 
accounted  for  the  most  variance  in  the  prediction.  Additionally,  the  analysis  led  to  construct 
models  by  munition  type.  This  resulted  in  four  models  as  shown  with  the  best  models,  1 A  and 
2A,  which  are  the  most  operationally  relevant. 
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[E3I  The  distribution  of  simulated  results  verses  predicted 
*32=0]  values  provide  for  risk  averse  adjustments _ 


Model  1 A  Model  IB  Model  2A  Model  2B 


a  -  Residual  sqrtPK  ^  ?  Residual  sqrtPK  2  A  -  Residual  sqrtPK  2 


a  -  Residual  sqrtPK 

01 


JFk  =  f(d.  r.  D0,Al,  Alo,  T.  0,  TLE ) 
JFK  =  nd.r,D0,AL,AL0.  T,  0,  TLE )  +  C 

Pk  =  [f(d,r,D0,AL,AL0lT,  0,TLE)  +  C]2 
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Model 

75%  Adj 

90%  Adj 

99%  Adj 

1A 

0.014 

0.032 

0.10 

IB 

0.04 

0.077 

0.15 

2A 

0.005 

0.011 

0.04 

2B 

0.022 

0.041 

0.10 

Significant  push  in  recent  conflicts  to  reduce  non-enemy  damage  (collateral  damage  and 
fratricide) 

In  order  to  meet  this  effort,  a  conservative  adjustment  can  be  made  to  better  protect  the 
safety  of  non-enemies 

Based  on  the  Model  in  use,  a  constant  based  on  error  distribution  can  be  added  to  the 
model  before  it  is  squared 

Example:  To  be  90%  certain  that  PK  is  this  value  or  less,  add  a  civilian  adjustment  C 
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[Ell  The  response  surface  modeling  provides  a  flexible 
approach  that  is  also  in  the  form  of  an  equation. 

Using  simulation  results  based  on  selected  factors,  using  the  Carleton 
Damage  Function,  quality  generalized  results  are  achieved. 


The  developed  model  provides  a  robust  and  flexible  capability  for  modeling 
effects  and  tradeoffs  between  precision  and  conventional  munitions. 


Methodology  is  adaptable  for  determining  killed,  wounded,  and  suppression 
effects. 


The  model  is  easily  modified  for  risk-averse  effects  on  non-combat 
populations. 
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[Ell  *n  summary>  both  models  provide  CAA  a  robust 
asea  precision-munitions  effects  modeling  capability. 

Our  research  team  developed  two  methods  to  model 
precision-guided  munitions  effects. 


The  analytical  models  improve  ultimately  on 
munitions  estimates  from  the  model  output. 


The  response  surface  models  provide  additional  risk 
analysis  capabilities  for  exploration. 


Questions? 
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Project  Team 


Sponsor:  US  Army  Munitions  Management  Division 
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Stakeholders: 
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[gj|  We  then  analyzed  optimal  aiming  points  of  M  weapons 

02=3]  in  the  presence  of  dependent  and  independent  errors. 

Optimal  aiming  points  for  maximizing  kill  probability: 
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We  analyzed  the  optimal  distribution  of  aim  points  initially  for  circular  and  elliptical 
shaped  aiming  areas.  The  series  of  graphics  on  this  slide  depict  our  final  analysis  of  aim 
position  distribution  without  constraint  to  the  target  area.  We  varied  the  number  of  weapons  to 
aim  at  a  target  located  at  (0,0)  from  a  minimum  of  five  weapons  to  a  maximum  of  18  weapon 
systems  firing  a  single  round  at  the  target. 
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[C3I  Also,  the  scaling  law  for  the  decay  of  survival 
SH  probability  showed  a  linear  fitting. _ 


log (— log ( 1 — Pun ))  VS.  logM 


is  approximately  linear. 
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To  find  a  phenomenological  fitting  to  the  decay  of  survival  probability  as  a  function  of  M, 
we  consider  the  form  of  1  —  pkm(M)  —  exp (— aM^).  If  the  survival  probability 
approximately  satisfies  this  relation,  then  the  plot  of  log[—  log  (1  —  pkiu )]  vs.  log (M)  would 
approximately  follow  a  straight  line.  In  other  words  log[—  log  (1  —  pkm)]  =  log(cr)  + 
/?log(M).  The  left  panel  shows  the  plot  is  very  close  to  a  straight  line. 

In  the  right  panel,  we  plot  pkm  vs.  M  and  the  fitting  function  1  —  exp(— 0.35 M0  73).  For 
the  set  of  parameters  values  used,  phenomenologically  we  have  the  approximation  of  pkiu  * 

1  -  exp(-0.35M°'73). 
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We  compare  the  results  of  optimization  over  constrained  patterns  with  those  of  the  overall 
optimization.  The  charts  show  plots  of  the  difference  in  optimal  (maximal)  damage  fractions 
(popt)  between  constrained  and  unconstrained  optimizations  as  a  function  of  area  target  radius 
(s)  for  various  number  of  weapons  (M).  The  difference  shows  how  accurate  it  is  to  optimize 
over  a  given  pattern.  For  M  <  6  (top  left  panel),  the  best  approximate  optimal  damage  fraction 
is  achieved  by  distributing  the  aiming  points  over  an  ellipse,  uniformly  in  polar  angle  (Pattern 
Bl).  At  M  =  7  (lower  left  panel),  the  best  approximate  optimal  damage  fraction  is  achieved  by 
placing  an  aiming  point  at  center  and  placing  the  rest  six  aiming  points  over  an  ellipse, 
uniformly  in  polar  angle  (Pattern  B2)  .  This  is  also  true  for  M  =  8  and  M  =  9.  As  the  number  of 
weapons  increases,  at  M  =  10  (top  right  panel),  the  best  approximate  optimal  damage  fraction 
is  achieved  by  placing  two  aiming  points  on  the  x-axis  and  the  rest  eight  aiming  points  over  an 
ellipse,  uniformly  in  polar  angle  (Pattern  B3).  The  constrained  optimum  over  Pattern  B3 
remains  very  accurate  at  M  =  12  weapons  (lower  right  panel). 
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To  explore  how  the  size  of  optimal  aiming  points  scales  with  the  area  target  radius,  we 
plotted  these  two  quantities  against  each  other  in  a  log-log  plot.  The  left  panel  compares  sets  of 
scaled  optimal  aiming  points  of  pattern  B3  for  s  =  100,  s  =  200,  s  =  300,  and  s  =  400.  The 
comparison  demonstrates  that  not  only  the  spread  of  optimal  aiming  points  scales  as  Vs  ,  the 
distribution  of  optimal  aiming  points  after  scaling  is  approximately  invariant  with  respect  to 
the  area  target  radius. 


The  right  panel  shows  that  the  size  of  the  optimal  aiming  points  is  approximately 
proportional  to  the  square  root  of  the  area  target  radius  (s).  This  observations  suggests  that  we 
should  normalize  the  aiming  points  by  the  square  root  of  the  area  target  radius. 
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Appendix  C 

Analytic  Solution  Methodology 


This  appendix  contains  the  technical  report  for  the  analytical  solution  of  bring  precision- 
guided  munitions  at  point  and  area  targets,  supporting  MATLAB  code  used  to  implement 
the  methodologies,  and  finally  a  white  paper  detailing  the  functions  produced  in  MATLAB. 
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Abstract 

We  study  the  damage  probability  when  M  weapons  are  used  against  a  unitary  tar¬ 
get.  We  use  the  Carleton  damage  function  to  model  the  distribution  of  damage  proba¬ 
bility  caused  by  each  weapon.  The  deviation  of  the  impact  point  from  the  aimpoint  is 
attributed  to  both  the  dependent  error  and  independent  errors.  The  dependent  error  is 
one  random  variable  affecting  M  weapons  the  same  way  while  independent  errors  are 
associated  with  individual  weapons  and  are  independent  of  each  other.  We  consider 
the  case  where  the  dependent  error  is  significant,  non- negligible  relative  to  independent 
errors.  We  first  derive  an  explicit  exact  solution  for  the  damage  probability  caused  by 
M  weapons  for  any  M.  Based  on  the  exact  solution,  we  find  the  optimal  aimpoint 
distribution  of  M  weapons  to  maximize  the  damage  probability  in  several  cases  where 
the  aimpoint  distribution  is  constrained  geometrically  with  a  few  free  parameters,  in¬ 
cluding  uniform  distributions  around  a  circle  or  around  an  ellipse.  Then,  we  perform 
unconstrained  optimization  to  obtain  the  overall  optimal  aimpoint  distribution  and  the 
overall  maximum  damage  probability,  which  is  carried  out  for  different  values  of  M,  up 
to  20  weapons.  Finally,  we  derive  a  phenomenological  approximate  expression  for  the 
damage  probability  vs  M,  the  number  of  weapons,  for  the  parameters  studied  here. 

Keywords:  Damage  probability;  Carleton  damage  function;  multiple  weapons  with  de¬ 
pendent  errors;  exact  solution;  optimal  distribution  of  aimpoint 


1  Introduction 

The  probability  of  killing  or  damaging  a  target  depends  heavily  on  how  close  a  weapon 
is  delivered  to  the  target.  This  delivery  accuracy  of  a  weapon  may  be  affected  by  many 
components.  In  general,  the  errors  are  usually  divided  into  two  main  groups:  the  dependent 
*  Corresponding  author,  hzhou@nps.edu 
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error  and  independent  errors.  The  dependent  error  is  related  to  the  aiming  error  that  results 
from  a  miscalculation  of  latitude,  longitude,  distance,  wind  effect,  or  uncertainty  in  locating 
the  target  position.  The  dependent  error  results  in  the  armament  impacting  away  from  the 
desired  target  point  and  it  affects  all  weapons  the  same  way.  The  independent  errors  refer 
to  ballistic  dispersion  errors,  which  may  result  from  variations  in  bullet  shape,  variations  in 
gun  barrels,  or  variations  in  amount  of  explosive  used  inside  each  bullet  [1]. 

Due  to  many  uncertainties  in  the  field  of  weapon  effectiveness,  Monte  Carlo  simulations 
have  been  widely  employed  to  estimate  the  probability  of  target  damage  [2].  Even  though 
Monte  Carlo  simulations  can  provide  reasonable  estimates,  exact  solutions  are  mathemat¬ 
ically  more  attractive  and  practically  more  useful.  The  objectives  of  this  paper  are  i)  to 
derive  explicit  exact  solution  for  the  damage  probability  caused  by  multiple  weapons  against 
a  single  target,  ii)  to  use  the  exact  solution  to  maximize  the  damage  probability  with  respect 
to  the  aimpoint  distribution  of  weapons,  with  or  without  geometric  constraint (s)  on  the 
aimpoint  distribution,  and  iii)  to  study  the  relation  of  damage  probability  to  the  number  of 
weapons  when  the  dependent  error  is  significant.  The  results  obtained  here  can  be  applied 
to  indirect  fire  artillery,  or  GPS/INS-guided  weapons. 

The  remainder  of  this  paper  will  progress  as  follows.  Section  II  provides  the  detailed 
mathematical  formulation  and  explicit  exact  solution  for  the  kill  probability.  Section  III 
considers  the  performances  of  various  aimpoint  distributions.  Finally,  Section  IV  presents 
conclusions  and  future  work. 


2  Mathematical  formulation 

We  consider  a  single  point  target  in  the  two  dimensional  space.  We  establish  the  coor¬ 
dinate  system  such  that  the  target  is  located  at  the  origin  point  Xtarget.  =  (0,0).  We  use  M 
weapons  with  dependent  and  independent  errors  to  fire  on  the  target.  Due  to  the  presence 
of  significant  dependent  error,  if  all  M  weapons  are  aimed  at  xtarget  =  (0,0),  the  M  impact 
points  may  be  uniformly  shifted  away  from  the  target  by  a  significant  distance,  resulting  in  a 
small  damage  probability.  To  make  the  damage  probability  less  susceptible  to  the  dependent 
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error,  we  aim  the  M  weapons  at  M  different  points  distributed  around  the  target.  When  the 
dependent  error  shifts  some  impact  points  away  from  the  target,  it  simultaneously  shifts  the 
some  other  impact  points  toward  the  target.  In  this  study  all  weapons  are  assumed  to  be 
perfectly  reliable.  Gross  errors  due  to  anomalies  such  as  catastrophic  weapon  system  failure, 
adverse  weapon  separation  effects,  and  GPS  jamming  are  neglected. 

Let 


•  fj  =  the  aiming  point  of  weapon  j. 

•  Y  =  miss  distance  from  the  aimpoint  due  to  the  dependent  error  of  M  weapons, 
affecting  the  impact  points  of  all  M  weapons  uniformly. 

•  Xj  =  miss  distance  from  the  aimpoint  due  to  the  independent  error  of  weapon  j, 
affecting  only  the  impact  point  of  weapon  j  individually.  We  assume  that  {X,-,  j  = 
1,  2, . . . ,  M}  are  independent  of  each  other  and  independent  of  random  variable  Y. 

The  impact  point  of  weapon  j  is  given  by 

Wj  =  rj  +  Y  +  Xj 

We  model  the  dependent  error  Y'as  a  normal  random  variable  with  zero  mean: 


where  o\  and  a 2  are  standard  deviations,  respectively,  in  the  two  coordinate  directions,  which 
give  an  indication  of  the  spread  of  the  dependent  error  in  the  two  directions.  We  model  each 
independent  error  Xj  as  a  normal  random  variable  with  zero  mean: 


Further,  we  assume  that  the  independent  errors  of  individual  weapons  {Xj,  j  =  1,  2, . . . ,  M} 
are  independent  of  each  other  and  are  independent  of  the  dependent  error  Y . 
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We  use  the  mathematical  fact  that  the  sum  of  two  independent  normal  random  variables 
is  a  normal  random  variable.  Suppose  U  r-u  iV(0,  a2)  and  V  ~  iV(0,  s2).  We  have 


Z  =  V  +  U  ~  iV(0,  s2  +  a2) 
The  probability  density  functions  of  1/  and  Id  are  given  by 


(1) 


Puiu) 

1 

—  pvn  i 

f-u2 

—  . -  U 

a/  27T  <t2 

\2a2 

pv(v ) 

1  1 

—  pvn  1 

f-v2' 

1 

k2s2  ' 

In  terms  of  the  probability  density  functions,  we  write  equation  (1)  as 


V2 


:  exp 


7TS 


-\z-u) 

2  s2 
1 


v/2 


:  exp 


7177- 


— M 
2^2 


du 


:  exp 


for  any  z 


y/M*2  +  *2)  \2(-S2  +  (72) 

Applying  a  change  of  variables  unew  =  —u,  denoting  unew  still  by  u  for  simplicity  and 
multiplying  the  equation  by  y/2ns2,  we  get 


exp 


-(z  +  uf 
2s 2 


a/27 T(72 


exp 


-u 


exp 


2a2 


du 


s2  +  a2  )  1  \2  (s2  +  a2) 

We  rewrite  the  equation  above  in  terms  of  expected  values: 


for  any  z 


E, 


u 


exp 


-(z  +  uy 

2s 2 


exp 


(2) 


s2  +  cr2  J  \2(s2  +  a2) 

Here  the  notation  Eu  indicates  the  average  with  respect  to  random  variable  U  while  z  and 
s2  are  fixed,  not  varying  with  U.  Equation  (2)  is  valid  for  any  normal  random  variable 
U  ~  iV(0,<72),  and  for  any  z  and  s 2  >  0.  In  the  analysis  below,  we  will  use  equation  (2) 
extensively. 

We  use  the  Carleton  damage  function  to  model  the  probability  of  killing  by  an  individual 
weapon.  Let  w  =  (w(l),  w(2))  be  the  impact  point  of  a  weapon  where  w(l)  and  w( 2)  describe 
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the  impact  points  in  the  range  and  deflection  directions  from  the  target.  The  probability  of 
the  target  being  killed  by  an  weapon  at  impact  point  w  is  modeled  mathematically  as 


Pr(target  being  killed  by  one  weapon  at  impact  point  w) 

= exp  exp  (3) 

This  is  called  the  Carleton  damage  function  or  the  diffuse  Gaussian  damage  function  [3]. 
The  two  parameters  b\  and  b2  hr  the  Carleton  damage  function  (3)  represent  the  effective 
weapon  radii  in  the  range  and  deflection  directions,  respectively.  With  the  impact  points  of 
the  M  weapons  given  by  Wj  =  fj  +  Y  +  Xj,  j  =  1,  2, ... ,  M,  the  probability  of  the  target 
located  at  the  origin  being  killed  by  the  M  weapons  is 


Pr(target  being  killed  by  M  weapons  at  impact  points  {wj,  j  =  1,  2, . . . ,  M}) 


all  k-combinations 


(.7 1 7  J  2 1  •  •  •  1  jk ) 
from  (1,2,...,  M) 


M 


e<-p  e  n«p 

k=1 


26? 


■  ,jk ) 


xneip(E!E|)tM) 

^ — - V - ' 

We  calculate  the  probability  of  the  target  being  killed  averaged  over  independent  errors 
{Xj,  j  =  1,  2, ... ,  N}  and  averaged  over  the  dependent  error  Y .  For  that  purpose,  we  only 
need  to  calculate  the  average  of  each  term  inside  the  summation:  E  [iy  (ji, . . .  ,jk)  x  F2(.7i,  •  •  • ,  jk)}- 
Notice  that  F\  (j { , . . .  ,jk)  involves  only  the  horizontal  components  and  F2 (j  1 , . . .  ,jk )  involves 
only  the  vertical  components  of  {Xj}  and  Y.  Since  the  horizontal  components  and  vertical 
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components  are  independent  of  each  other,  we  have 


E  . . .  ,jk )  x  F2(jh  . . .  ,jk)]  =  E  . . .  ,jk)]  x  E  [F2(j i, . . .  ,jk)] 

Since  F\  (j i , . . .  ,jk)  and  F2 (j  1 , . . .  ,jk)  have  exactly  the  same  format,  we  only  need  to  derive 
the  analytical  expression  for  one.  For  conciseness,  we  denote  7^(1),  Y(  1),  and  Xjt(  1)  simply 
by  r jt ,  Y,  and  Xjt  in  the  calculation  of  E  [Fi(ji, . . .  ,jk)].  We  first  average  F\  (j j , . . . ,  jk)  over 
independent  errors  {Xj  rs_/  N(0,dl)}. 


E{Xj}  [Fi(ji,  ■  ■  ■  ,jk )]  -  E{Xj}  n  exp 


~(rji + x. + 
2  6? 


/rn  r  /-(r?f  +Y  +  XV)2 

fl  E*n 


Each  term  in  the  product  is  an  average  of  the  form  on  the  left  hand  side  of  (2).  Applying 
equation  (2),  we  write  each  average  as 


Exv  exp 


~{r 7,  +Y-\-  Xj.)2 


bj  y  f-Y+Y? 

W^J  exp[~W[YY) 


Substituting  this  result  into  equation  (5),  we  obtain 

E ™  |F'fe'  -  -  -  'i)] = n  (ifTif) 2  exp  (yy 


6?  +  d? 


EliY  +  Yf_ 
2(6?  +  d\) 


bj  V 
bj  +  dj) 


-(Y*  +  2YYL1rjl/k  +  YL1rl/k 

2(6?  +  dj)/k 


WTdi)  exp  1 - 2(6?  +  d?) /k -  exp 


-  [Y  +  £f=i  r3i/k 
2(6?  +  dj)/k 


Next  we  average  over  the  dependent  error  Y  ~  N(0,af).  Again,  the  average  is  of  the  form 
on  the  left  hand  side  of  (2).  Applying  equation  (2),  we  arrive  at 


EY  exp 


(y  +  Eti  rjjk)  \  /  (6 ?  +  rf?)/fc  y  I  -(Eti^A) 

2(bj  +  dj)/k  \(6?  +  dj)/k  +  <j\)  1  2((6?  +  d?)//c  +  a\ 
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Thus,  the  overall  average  of  Fx  (j‘ x , . . . ,  jk)  has  the  expression 


bl 


(bi  +  di)/k  y 


x  exp 


b\  +  d\)  \(bl  +  df)/k  +  af 


2(6?  +  d\)/k 


ml  +  <*?)/*  +  < 


Similarly,  the  overall  average  of  F2  (j x , . . .  ,jj  has  the  expression 

{bl  +  df)/k 


E  [F2(j'i,  . . .  ,jk)\  = 

/ 

V 


bj  +  dl 


+  d%)/k  + 


x  exp 


Ef.i  n(2)/*)2  -  Et,  o,(2)2/fc  (eLi  o,(2)/fc); 


\ 


2(6?  +  d|)/fc 


2((&2  +  +  a2 ) 


J 


(6) 


(7) 


The  probability  of  target  being  killed,  averaged  over  independent  errors  and  dependent  error, 
is  called  kill  probability,  and  is  denoted  by  P^[\\{M  weapons).  It  has  the  expression 

M 

Pkill(M  weapons)  =  -  J^(-l)fc  E  [F^Uu  •  •  •  Jk)]  E  [F2(ju  ■  ■  ■  Jk)]  (8) 

fc=1  (iii-*) 


where  E  [Fjji, . . . ,  jj]  and  E  [F2(ji, . . . ,  jj]  are  given  in  (6)  and  (7)  above.  Together, 
equations  (6),  (7),  and  (8),  give  us  an  explicit  analytical  expression  for  calculating  the  kill 
probability. 


After  the  completion  of  the  above  derivation,  we  discovered  that  similar  approaches  had 
been  taken  separately  by  von  Neumann  [4]  and  by  Washburn  [5]. 


3  Performances  of  various  aimpoint  distributions  of 
multiple  weapons  against  a  single  target 


Now  we  apply  the  exact  solution  to  examine  the  kill  probability  corresponding  to  various 
distributions  of  the  aimpoints  of  M  weapons. 
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Let  Al  denote  the  weapon  lethal  area  or  the  fragmentation  mean  area  of  the  effectiveness. 
It  describes  the  effect  of  a  warhead  against  a  target  and  includes  the  effects  of  direct  hit, 
blast,  and  fragmentation.  We  can  calculate  Al  from  the  Carleton  damage  function  (3)  as 


f*  oo  poo 


AL=LJ-Jxp{~i-i]dxdy=Mib2 


(9) 


The  aspect  ratio  of  the  weapon  radii  of  the  Carleton  damage  function  a  =  ^  is  described  by 
the  empirical  formula: 

a  =  max(l  —  0.8  cos  9,  0.3)  (10) 


where  9  is  the  impact  angle. 


Once  the  lethal  area  Al  and  the  aspect  ratio  a  are  given,  one  can  calculate  the  weapon 
radii  for  the  Carleton  damage  function  (3)  as  follows: 


bi  = 

b2  = 


(11) 

(12) 


For  all  the  cases  considered  in  this  paper,  we  choose  Al  =  2270  ft2,  9  —  65°.  This  yields 
bi  =  15.4640  ft  and  b2  =  23.3628  ft.  Furthermore,  we  choose  g\  —  g2  —  30  for  the  dependent 
error  and  d\  =  d2  =  5  for  the  independent  errors. 


We  first  consider  the  case  of  M  weapons  with  aimpoints  uniformly  distributed  on  a  circle 
as  formulated  below 


fj  =  r  cos  9  + 


2?r (j  ~  1) 
M 


,  r  sin  d  + 


2?r(j  ~  1) 
M 


where  r  is  the  radius  and  9  the  phase  off-set  angle  of  the  distribution.  These  are  parameters 
that  we  can  tune  to  maximize  the  kill  probability. 


For  each  value  of  M,  we  maximize  the  kill  probability  with  respect  to  (r,  9).  This  un¬ 
constrained  nonlinear  optimization  can  be  achieved  by  using  MATLAB  built-in  function 
fminsearch  which  is  based  on  a  direct  search  method  of  Lagarias  et  al.  [6].  The  results  are 
listed  in  Table  1. 
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Table  1:  The  optimal  distribution  for  M  aimpoints  when  they  are  uniformly  distributed 
around  a  circle  and  the  corresponding  probability  of  kill.  Here  r  is  the  radius  and  6  is  the 
phase  off-set  angle. 


M 

ropt 

^opt 

Ckill 

1 

0 

*** 

0.27597 

2 

16.246 

0 

0.43690 

3 

22.960 

7T 

6 

0.53834 

4 

26.948 

0 

0.62291 

5 

29.192 

7 T 

To 

0.68212 

6 

31.086 

7 r 

6 

0.72869 

7 

32.529 

7 r 

14 

0.76474 

8 

33.731 

0 

0.79360 

9 

34.747 

7 r 

18 

0.81702 

10 

35.63 

7 T 

To 

0.83635 

11 

36.409 

7T 

22 

0.85251 

12 

37.105 

0 

0.86617 

Note  that  the  Carleton  damage  function  we  use  is  not  isotropic.  It  has  different  effective 
radii  in  the  range  and  deflection  directions.  To  accommodate  this  anisotropic  property  of  the 
Carleton  damage  function,  we  consider  the  case  of  M  weapons  with  aimpoints  distributed 
on  an  ellipse  as  formulated  below 


ri  = 


(qVv 


cos 


+ 


2vr(j-H)\  _q_ 

M  )'  ,/r} 


Sill 


+ 


Mj  - 1) 

M 


where  r)  is  the  aspect  ratio  of  the  ellipse.  In  the  formulation  above,  we  elongate  one  axis  by 
y/rj  and  simultaneously  shrink  the  other  axis  by  the  same  factor.  In  this  way,  the  area  of  the 
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Table  2:  The  optimal  distribution  for  M  aimpoints  when  they  are  uniformly  distributed 
around  an  ellipse  and  the  corresponding  probability  of  kill.  Here  0  is  the  off-set  value  in 
the  parametric  equation  of  the  ellipse.  The  cases  of  M  —  1  and  M  =  2  are  not  affected  by 
aspect  ratio. 


M 

(major  axis)opt 

(minor  axis)Qp^ 

^opt 

hkill 

1 

0 

0 

*** 

0.27597 

2 

16.246 

16.246 

0 

0.43690 

3 

25.637 

17.639 

7 r 

6 

0.53989 

4 

29.621 

23.068 

0 

0.62477 

5 

30.411 

27.235 

7 r 

To 

0.68264 

6 

32.859 

28.548 

7 r 

6 

0.72958 

7 

34.292 

30.095 

7T 

14 

0.76560 

8 

35.848 

30.967 

0 

0.79469 

9 

37.135 

31.75 

7T 

18 

0.81829 

10 

38.342 

32.349 

7T 

To 

0.83784 

11 

39.436 

32.861 

7T 

22 

0.85420 

12 

40.457 

33.29 

0 

0.86806 

ellipse  is  maintained  at  n g2,  independent  of  the  aspect  ratio  77.  Parameter  q  has  the  meaning 


'area  of  ellipse  r - r - - - r 

q  =  \j - =  y  (major  axis)  x  (minor  axis) 


7T 


From  q  and  77,  we  can  determine  the  major  axis  and  the  minor  axis  as 


major  axis  =  qy/rj 

q 

minor  axis  =  - 

Vv 

We  should  point  out  that  parameter  (j)  is  not  the  polar  angle  of  the  aimpoint  of  weapon  1.  </> 
is  the  angular  value  used  in  the  parametric  equation  of  the  ellipse  to  calculate  the  aimpoint 
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of  weapon  1.  0  is  the  phase  angle  before  the  major  axis  is  elongated  and  before  the  minor 
axis  is  shrunk. 

Table  3:  The  optimal  distribution  for  M  aimpoints  when  one  of  them  is  aimed  at  the  origin 
while  the  rest  of  aimpoints  are  uniformly  distributed  around  an  ellipse,  and  the  corresponding 
probability  of  kill.  Here  0  is  the  off-set  value  in  the  parametric  equation  of  the  ellipse.  For 
the  cases  of  M  <  5,  the  kill  probability  is  not  improved  by  moving  one  of  the  M  aimpoints 
to  the  center. 


M 

(major  axis)opt 

(minor  axis)0p^ 

‘/’opt 

dkill 

1 

*** 

*** 

*** 

0.27597 

2 

22.161 

22.161 

0 

0.40957 

3 

25.412 

25.412 

0 

0.53737 

4 

32.918 

23.814 

7T 

6 

0.60947 

5 

34.369 

30.581 

0.1407vr 

0.67798 

6 

36.213 

33.451 

7T 

To 

0.73052 

7 

38.374 

34.765 

0 

0.77123 

8 

39.45 

36.17 

7T 

14 

0.80221 

9 

40.859 

36.86 

7T 

8 

0.8274 

10 

41.814 

37.655 

71 

18 

0.84766 

11 

42.838 

38.163 

0 

0.86449 

12 

43.709 

38.648 

71 

22 

0.87853 

For  each  value  of  M,  we  maximize  the  kill  probability  with  respect  to  (q,  0, 77).  We  obtain 
the  results  in  Table  2. 

In  the  above,  we  calculated  the  performance  of  placing  the  aimpoints  of  M  weapons  along 
a  circle  or  an  ellipse.  We  now  examine  the  case  of  aiming  one  weapon  at  the  center  and 
aiming  the  rest  (M  —  1)  weapons  at  positions  distributed  on  an  ellipse.  The  aimpoints  of  M 
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weapons  are  distributed  as  formulated  below. 


r-i  =  I  qy/rj  cos  |  0  + 


j 


M  —  1  )  y/fj  V  M  ~  1 


,  3  =  1,2, 


M-  1 


fM  =  (0, 0) 


For  each  value  of  M,  we  maximize  the  kill  probability  with  respect  to  (q,  0,  rj).  The 
optimal  results  are  reported  in  Table  3. 

Next,  we  fully  optimize  the  distribution  of  M  aimpoints  without  constraining  them  on  a 
circle  or  an  ellipse.  We  represent  the  M  aimpoints  in  polar  coordinates. 


(O,  9j)  ,  j  =  l,2,...,M 

The  optimal  solutions  for  M  —  1,  M  =  2,  M  =  3  and  M  —  4  are  listed  in  Table  4. 

Figure  1  shows  the  optimal  distributions  of  aimpoints  for  M  —  1  (yellow  circles)  and 
M  =  2  (blue  squares)  while  the  optimal  distributions  for  the  cases  of  M  =  3  (yellow  circles) 
and  M  =  4  (blue  squares)  are  displayed  in  Figure  2. 


O  M  =  1 
■  M  =  2 


Figure  1:  Optimal  distributions  of  aimpoints  for  M  =  1  (yellow  circles)  and  M  =  2  (blue 
squares) . 
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O  M  =  3 
■  M  =  4 


II 


Figure  2:  Optimal  distributions  of  aimpoints  for  M  =  3  (yellow  circles)  and  M  =  4  (blue 
squares) . 


O  M  =  5 
■  M  =  6 


Figure  3:  Optimal  distributions  of  aimpoints  for  M  =  5  (yellow  circles)  and  M  =  6  (blue 
squares) . 


The  optimal  solutions  for  M  —  5,  M  —  6,  M  —  7  and  M  =  8  are  given  in  Table  5. 
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Table  4:  Optimal  distributions  of  aimpoints  and  the  corresponding  probabilities  of  kill  for 
M  =  1,  M  =  2,  M  =  3  and  M  =  4. 


M  =  1 

M  =  2 

M  =  3 

M  =  4 

Pkiu  =  0.27597 

Pkiii  =  0.4369 

Pkiii  =  0.53989 

Pkiii  =  0.62477 

3 

ri 

ri 

h 

r3 

n 

0j 

1 

0 

16.246 

0 

23.975 

0.122117T 

29.621 

0 

2 

16.246 

7T 

23.975 

0.87789tt 

23.068 

0.57T 

3 

17.411 

1.57T 

29.621 

71 

4 

23.068 

2vr 

Table  5:  Optimal  distributions  of  aimpoints  and  the  corresponding  probabilities  of  kill  for 
M  =  5,  M  =  6,  M  =  7  and  M  =  8. 


M  =  5 

M  =  6 

M  =  7 

M  =  8 

Pkiii  =  0.68505 

Pkiu  =  0.73391 

Pkiii  =  0.77218 

Pkiii  =  0.80341 

3 

r3 

ri 

*3 

r3 

*3 

r3 

*3 

1 

33.839 

0.08879tt 

38.334 

0 

40.69 

0 

43.43 

0 

2 

24.734 

0.57T 

27.589 

0.33554vr 

34.456 

0.32695tt 

36.016 

0.28566vr 

3 

33.839 

0.911217T 

27.589 

0.664467T 

34.456 

0.67305tt 

34.191 

0.5809tt 

4 

26.353 

1.3057tt 

38.334 

71 

40.69 

71 

40.507 

0.85228vr 

5 

26.353 

1.6943tt 

27.589 

1.3355tt 

34.456 

1.327tt 

40.507 

1.1477tt 

6 

27.589 

1.6645tt 

34.456 

1.673tt 

34.191 

1.41917T 

7 

0 

36.016 

1.7143tt 

8 

0.6587 

0 

Figure  3  illustrates  the  optimal  distributions  of  aimpoints  for  M  =  5  (yellow  circles)  and 
M  —  6  (bine  squares);  Figure  4  shows  the  optimal  distributions  of  aimpoints  for  M  —  7 
(yellow  circles)  and  M  —  8  (bine  squares). 

The  optimal  solutions  for  M  =  9 ,  M  —  10,  M  —  11  and  M  =  12  are  listed  in  Table  6. 
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Figure  4:  Optimal  distributions  of  aimpoints  for  M  =  7  (yellow  circles)  and  M  —  8  (blue 
squares) . 


Figure  5:  Optimal  distributions  of  aimpoints  for  M  =  9  (yellow  circles)  and  M  =  10  (blue 
squares) . 

Figure  5  displays  the  optimal  distributions  of  aimpoints  for  M  —  9  (yellow  circles)  and 
M  —  10  (blue  squares);  the  optimal  distributions  for  M  =  11  (yellow  circles)  and  M  =  12 
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Table  6:  Optimal  distributions  of  aimpoints  and  the  corresponding  probabilities  of  kill  for 
M  =  9,  M  =  10,  M  =  11  and  M  =  12. 


M  =  9 

M  =  10 

M  =  11 

M  =  12 

Pkui  =  0.82935 

Pkui  =  0.85158 

Pkui  =  0.86957 

Pkui  =  0.88499 

3 

rj 

% 

rj 

ri 

rj 

1 

46.336 

0 

48.534 

0 

50.078 

0.05606?r 

49.091 

0.189727T 

2 

36.994 

0.24837vr 

42.116 

0.25814?r 

42.093 

0.28098?r 

41.676 

0. 40044 7T 

3 

35.587 

0.57T 

38.162 

0.57T 

38.939 

0.57T 

41.675 

0.59952?r 

4 

36.994 

0.75163vr 

42.116 

0.741867T 

42.093 

0.71902? r 

49.089 

0.81023?r 

5 

46.336 

71 

48.534 

71 

50.078 

0.943947T 

52.524 

1.062?r 

6 

36.994 

1.2484? r 

42.116 

1.2581?r 

45.328 

1.1852?r 

45.409 

1.296?r 

7 

35.587 

1.57T 

38.162 

1.57 r 

38.882 

1.395?r 

40.554 

1.57 r 

8 

36.994 

1.7516? r 

42.116 

1.7419?r 

38.882 

1.605?r 

45.407 

1.93  797T 

9 

0 

12.452 

0 

45.328 

1.8148?r 

52.523 

1.7578?r 

10 

12.452 

71 

12.535 

0.023727T 

21.097 

1.9977?r 

11 

12.535 

0.976287T 

0.44808 

1.5002?r 

12 

21.097 

1.0023?r 

(bine  squares)  are  plotted  in  Figure  6. 

As  M  (the  number  of  weapons)  increases,  the  optimal  distribution  of  aimpoints  has  more 
layers,  covering  a  larger  area  with  a  more  uniform  distribution  over  the  area.  In  Figure  7, 
we  plot  the  optimal  distributions  of  aimpoints  for  M  =  15  (yellow  circles)  and  M  —  18  (blue 
squares) . 

Next,  we  study  the  optimal  kill  probability  as  a  function  of  M.  Let  Pkui(M)  denote  the 
kill  probability  corresponding  to  the  optimal  distribution  of  aimpoints  for  the  case  of  M 
weapons.  As  M  increases,  the  survival  probability  of  the  target,  (1  — Pkui(M )),  decreases. 
In  the  absence  of  dependent  error  and  when  the  aimpoints  are  all  fixed  at  one  point,  the 
outcome  of  each  weapon  affected  by  its  independent  error  is  statistically  independent  of  the 
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Figure  6:  Optimal  distributions  of  aimpoints  for  M  =  11  (yellow  circles)  and  M  —  12  (blue 
squares) . 


Figure  7:  Optimal  distributions  of  aimpoints  for  M  —  15  (yellow  circles)  and  M  —  18  (blue 
squares) . 

outcome  of  other  weapons  affected  by  their  own  independent  errors.  In  this  situation,  the 
probability  of  surviving  M  weapons  is  simply  the  M- th  power  of  the  probability  of  surviving 
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one  weapon:  1  —  Pkm(M)  =  (1  —  Pkiii(^))M ■  In  other  words,  in  the  absence  of  dependent 
error,  the  log  survival  probability  is  a  linear  function  of  M. 

log(l  -pkui(M ))  =  M  log(l  -pkui(  1)) 


In  the  presence  of  dependent  error,  however,  the  situation  is  completely  different.  The  same 
dependent  error  affects  all  M  weapons.  The  outcomes  of  individual  weapons  are  no  longer 
independent  of  each  other.  As  a  matter  of  fact,  when  the  M  weapons  are  all  aimed  at  the 
same  position,  the  outcomes  of  individual  weapons  are  highly  correlated  with  each  other.  As 
an  example,  we  examine  the  case  of  aiming  all  M  weapons  at  the  origin.  The  averages  of 
. . .  ,jk)  and  F2(j i, . . . ,  jk)  are  calculated  from  equations  (6)  and  (7)  as 


E[F1(J1,...,jk)] 


bj  V  (  ( bl  +  dp/h  y 
b\  +  d\)  \(bl  +  df)/k  +  o\  ) 


E[F2(j1,...,jkj) 


(  M  \ 1  (  (Pl  +  dp/k  \  * 
\b22  +  dl)  \(bl  +  d22)/k  +  a\) 


The  kill  probability  is 


M 

Pmi{M)  =  ~Yy  1)"  E  E\F1(jll...,jk)\E\F2(j1,...,jk)\ 
k=1 


M 


k= 1 


M 


bib-2 


b\  +  d\ 


bl  +  dl 


k  /  \  K h 2  772Y 7p  1 

h '  /  V  V  +  «iK°2  d 


+  df  +  ka'f  V  +  fco-l 


In  the  absence  of  dependent  error,  we  have  a\  —  a2  —  0,  and  the  kill  probability  is 

\  M 

1  -  Pkui(M)  =  [  1  blb 2 


V  (bj  +  d\){b2  +  d|)i 

In  the  presence  of  dependent  error,  to  simplify  the  analysis,  we  assume  that  the  independent 


cr?  (To 


errors  are  zero  (<A  =  d2  =  0)  and  assume  that  — f  =  -A  =  ca2.  The  kill  probability  becomes 

b\  b2 


M 


l-pkill{M)  =  1  + 

For  the  first  few  values  of  M,  we  obtain 


k= 1 


M 


k  J  1  +  /eta2 


1  —  Pfcizz(l)  = 


ta 


1  +ta2 
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1  —  Pkill( 2) 


f  u2  \  (  2oo2 
Vl  +  wV  \  1  +  2a;2 

f  uj2  \  f  2cn2 
Pkiui  )  -  ^1  +  w2J  ^1  +  2^2 

Using  mathematical  induction,  we  can  prove  that 

i-p““(M)  =  n(rrb) 

k= i  v  7 

Clearly,  when  all  M  weapons  are  aimed  at  the  same  positon,  1  —  Pkm(M)  decays  less  than 
geometrically  with  M. 


T3 

0 


6 


a, 


M  (number  of  weapons) 


Figure  8:  Left  panel:  Comparison  in  the  decay  of  survival  probability  (1  —pkui),  of  the  case  of 
aiming  all  M  weapons  at  (0,  0)  vs  the  case  of  using  optimal  distribution  of  aimpoints.  Right 
panel:  Enhancement  in  the  decay  of  survival  probability  (1  —  Pkui)  attributed  to  optimizing 
the  distribution  of  M  aimpoints. 

With  the  optimal  distribution  of  aimpoints  for  M  weapons,  we  may  expect  that  1  — 
Pkui(M)  decays  faster  than  in  the  case  of  aiming  all  M  weapons  at  the  same  position.  Indeed, 
as  demonstrated  in  the  left  panel  of  Figure  8,  when  the  M  weapons  are  aimed  according 
to  the  optimal  distribution  of  aimpoints,  1  —  Pkui(M)  decays  much  faster  than  in  the  case 
of  aiming  all  M  weapons  at  (0,0).  The  right  panel  of  Figure  8  shows  the  enhancement 
in  the  decay  of  survival  probability  (1  —  Pkui(M))  attributed  to  the  optimal  distribution 
of  aimpoints.  Specifically,  in  the  right  panel  of  Figure  8,  we  plot  the  quantity  below  as  a 
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Figure  9:  Left  panel:  plot  of  log(l  —pkiii)  vs  M .  Right  panel:  plot  of  log(l  —pkui)  vs  log (M) 


function  of  M 


1  {all  aimed  at  (0,0)} 

1  ~~  Pkill 


-%  {optimal  aim  positions} 

1  —  Pkill 


M  (number  of  weapons) 


Figure  10:  Left  panel:  plot  of  log[—  log(l  —  Pkui)]  vs  log(M).  Right  panel:  plot  of  pkm  vs  M 
Even  with  the  optimal  distribution  of  aimpoints,  however,  the  log  survival  probability, 
log(l  —  Pkui(M)),  does  not  decrease  linearly  with  respect  to  M  in  the  presence  of  dependent 
error.  In  the  left  panel  of  Figure  9,  we  plot  log(l  —  pkui)  vs  M.  It  is  clear  that  in  the  presence 
of  dependent  error,  the  survival  probability  decreases  slower  than  the  geometric  decay. 

After  excluding  the  geometric  decay,  we  explore  the  possibility  of  a  power  law  decay  for 
the  survival  probability.  Specifically  we  examine  whether  or  not  the  survival  probability 
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obeys  the  power  law  1  —pkui(M)  =  aM  13 .  If  the  survival  probability  follows  this  power  law 
relation,  then  the  plot  of  log(l  —  Pkui)  vs  log(M)  would  be  a  linear  function 

log(l  -  pkui)  =  log  (a)  -  log  (M) 

In  the  right  panel  of  Figure  9,  we  plot  log(l  —pkui)  vs  log(M).  The  plot  demonstrates  clearly 
that  the  survival  probability  does  not  follow  a  power  law  decay. 

To  find  a  phenomenological  fitting  to  the  decay  of  survival  probability  as  a  function  of  M, 
we  consider  the  form  of  1  —pkui(M)  =  exp(— o :M/3).  If  the  survival  probability  approximately 
satisfies  this  relation,  then  the  plot  of  log[—  log(l  —  Pkui)]  vs  log(M)  would  approximately 
follow  a  straight  line. 


log[-  log(l  -  Pkui )]  =  log  (a)  +  p  log  (M) 

In  the  left  panel  of  Figure  10,  we  plot  log[—  log(l  —  pkm)\  vs  log(M).  The  plot  is  very  close 
to  a  straight  line.  In  the  right  panel  of  Figure  10,  we  plot  Pkm  vs  M  and  the  fitting  function 
1  —  exp(— 0.35M073).  For  the  set  of  parameter  values  used,  phenomenologically  we  have  the 
approximation: 

Pkui  «1-  exp(— 0.35M0'73) 


4  Conclusion 

We  have  considered  the  damage  probability  caused  by  multiple  weapons  against  a  single 
target.  Explicit  exact  solution  was  derived  for  the  damage  probability  in  the  case  of  M 
weapons  with  both  dependent  error  and  independent  errors.  Then  we  applied  the  explicit 
exact  solution  to  maximize  the  damage  probability  and  find  the  corresponding  optimal  dis¬ 
tribution  of  aimpoints.  We  observed  that  in  the  presence  of  significant  dependent  error,  the 
decay  of  the  survival  probablty  corresponding  to  the  optimal  aimpoints  distribution  (i.e. ,  1 
—  optimal  damage  probability)  is  slower  than  the  exponential  decay  with  respect  to  M,  the 
number  of  weapons.  This  observation  demonstrates  that  increasing  M  is  much  less  effective 
in  overcoming  the  dependent  error  than  in  overcoming  independent  errors.  We  found  that 
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phenomenologically  the  survival  probability  decays  exponentially  with  respect  to  a  fractional 
power  of  M.  Presumably,  the  fraction  power  varies  with  the  parameter  values  of  the  prob¬ 
lem.  The  mathematics  behind  this  phenomenological  expression  and  the  dependence  of  the 
fraction  power  on  the  parameter  values  will  be  investigated  in  future  studies. 


5  Disclaimer 

H.  Zhou  would  like  to  thank  TRAC-M  for  supporting  this  work.  The  views  expressed  in 
this  document  are  those  of  the  authors  and  do  not  reflect  the  official  policy  or  position  of 
the  Department  of  Defense  or  the  U.S.  Government. 
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Abstract 

This  paper  investigates  the  effect  of  launching  multiple  weapons  against  an  area 
target  of  normally  distributed  elements.  We  provide  an  analytical  form  of  the  average 
damage  fraction  and  then  apply  it  to  obtain  optimal  aimpoints.  To  facilitate  the 
computational  efforts  in  practice,  we  also  consider  optimizations  over  given  constrained 
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patterns  of  aimpoints.  Finally,  we  derive  scaling  laws  for  optimal  aimpoints  and  optimal 
damage  fraction  wtih  respect  to  the  radius  of  the  area  target. 

Military  OR  application  area:  Strike  Warfare 

OR  methodology:  Probabilistic  Operations  Research;  Decision  in  the  presence  of  uncer¬ 
tainty 


1  Introduction 

The  theory  of  firing,  which  mainly  concerns  about  aiming,  kill  probability  and  allocation 
of  munitions,  was  inspired  by  World  War  II  and  has  been  progressed  significantly  in  the  past 
decades  [1] .  A  brief  history  of  firing  theory  can  be  found  in  Washburn  and  Kress’s  book  [2] 
where  the  authors  also  presented  a  detailed  discussion  on  shooting  without  feedback  or  with 
feedback.  Another  good  reference  on  weaponeering  is  given  by  Driels  [3]. 

In  this  paper  we  are  interested  in  studying  the  effect  of  precision-guided  munitions  such  as 
Excaliburs.  These  coordinate-seeking  munitions  are  usually  guided  by  radio,  radar,  or  laser 
and  launched  by  a  cannon.  They  are  intended  to  hit  a  target  accurately  and  cause  minimal 
collateral  damage  to  civilians,  friendly  forces  and  infrastructure,  especially  hospitals,  schools, 
churches,  and  residential  homes.  The  precision-guided  weapons  are  in  general  subject  to 
target-location  errors  and  ballistic  dispersion  errors.  The  target-location  errors,  or  aiming 
errors,  result  from  inaccuracies  associated  with  identifying  a  target’s  location.  In  contrast, 
the  ballistic  dispersion  errors  are  caused  by  random  weapons  effects,  which  may  vary  from 
one  weapon  to  another  and  are  assumed  to  be  independent  from  shot  to  shot.  When  a  single 
weapon  is  fired,  it  is  natural  to  aim  it  at  the  expected  center  of  the  target.  However,  when 
multiple  weapons  are  launched  against  a  unitary  target,  the  probability  of  damaging  the 
target  can  be  improved  significantly  by  spreading  the  aimpoints  around  the  target  and  the 
optimal  distribution  of  aimpoints  has  been  investigated  in  our  recent  work  [4],  Our  goal  here 
is  to  extend  our  previous  studies  to  estimate  the  probability  of  destroying  an  area  target  of 
normally  distributed  elements  with  multiple  weapons.  We  will  seek  optimal  aimpoints  for 
various  number  of  weapons. 
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The  plan  of  this  paper  is  to  first  review  our  previous  analytical  results  for  the  case  of 
multiple  weapons  against  a  single  target  in  section  2.  Section  3  introduces  the  mathematical 
problem  of  multiple  weapons  being  released  against  an  area  target  consisting  of  normally 
distributed  elements.  Exact  solution  for  the  average  damage  fraction  is  then  derived.  Section 
4  calculates  the  optimal  aiming  points  and  examines  the  relation  among  the  radius  of  area 
target,  the  number  of  weapons  and  the  optimal  (maximum)  damage  fraction.  In  addition  to 
the  unconstrained  overall  optimization  of  the  damage  fraction,  we  also  study  empirical,  fast 
and  robust  constrained  optimization  over  several  prescribed  patterns.  The  goal  is  to  reduce 
the  computational  complexity  of  optimization  and  to  compute  a  set  of  nearly  optimal  aiming 
points  efficiently.  Section  5  provides  scaling  laws  for  optimal  aimpoints  and  optimal  damage 
fraction  with  respect  to  the  radius  of  area  target.  Section  6  highlights  conclusions. 


2  Review  of  our  previous  analytical  results  for  the  case 
of  multiple  weapons  against  a  single  target 

Even  though  the  world  is  three-dimensional,  most  targets  are  known  to  be  on  the  surface 
of  the  Earth  and  therefore  the  targets  are  assumed  to  be  in  a  two-dimensional  ground  space. 
Conventionally,  we  use  two  coordinates  to  define  this  ground  plane:  the  range  direction  and 
the  deflection  direction.  The  range  direction  is  defined  by  the  direction  of  the  weapon’s 
velocity  vector,  whereas  the  deflection  direction  is  perpendicular  to  the  range  direction. 

Previously  [4]  we  have  studied  the  case  of  multiple  weapons  with  both  dependent  and 
independent  errors  against  a  single  target  at  xtarget  =  (0,0).  For  reader’s  convenience,  we 
review  briefly  the  mathematical  formulation  of  the  problem.  Let 

•  fj  =  the  aiming  point  of  weapon  j  for  j  =  1,  2, ...,  M. 

•  Y  =  the  dependent  error  of  M  weapons,  affecting  the  impact  points  of  all  M  weapons 
uniformly.  We  assume  Y  is  a  normal  random  variable. 

•  Xj  =  independent  error  of  weapon  j,  affecting  only  the  impact  point  of  weapon  j 
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individually.  We  assume  that  {Xj,  j  =  are  normal  random  variables, 

independent  of  each  other  and  independent  of  normal  random  variable  Y. 


We  model  the  dependent  error  Y'as  a  normal  random  variable  with  zero  mean: 


where  0\  and  a2  are  standard  deviations,  respectively,  in  the  range  and  the  deflection  direc¬ 
tions,  which  give  an  indication  of  the  spread  of  the  dependent  error  in  these  two  directions. 
We  model  each  independent  error  Xj  as  a  normal  random  variable  with  zero  mean: 


The  impact  point  of  weapon  j  is  given  by 

Wj  =  fj  +  Y  +  Xj 

We  use  the  Carleton  damage  function  to  model  the  probability  of  killing  by  an  individual 
weapon.  Let  w  =  (w(l),  w(2))  be  the  impact  point  of  a  weapon  where  w(l)  and  w( 2)  describe 
the  impact  point  in  the  range  and  deflection  directions  from  the  target.  The  probability  of 
the  target  being  killed  by  an  weapon  at  impact  point  w  is  modeled  mathematically  as 


Pr(target  being  killed  by  one  weapon  at  impact  point  w) 

=exp(=1^)exp(=i|L)  (1) 

This  is  the  well-known  Carleton  damage  function  or  the  diffuse  Gaussian  damage  function 
[2],  The  two  parameters  b\  and  b2  in  the  Carleton  damage  function  (1)  represent  the  effective 
weapon  radii  in  the  range  and  deflection  directions,  respectively. 

The  probability  of  a  target  being  killed  by  the  M  weapons,  averaged  over  all  random  errors 
(i.e.,  dependent  and  independent  errors),  is  called  the  kill  probability  and  is  mathematically 
denoted  by  (target,  M  weapons).  Note  that  in  the  notation  for  the  kill  probability,  the 
target  identity  is  explicitly  included.  This  will  be  very  convenient  later  in  the  discussion  of 
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an  area  target  with  multiple  target  elements,  in  which  we  can  study  the  kill  probability  for 
each  individual  element. 

With  the  impact  points  of  the  M  weapons  given  by  {wj  —  fj  +Y  +  Xj,  j  =  1,2 , ,M}, 
we  derived  an  analytical  expression  for  the  kill  probability  as  a  function  of  quantities 
(01,  ^2,  di,  d2,  bi,  b2,  {fj,  j  —  1,2, ,  M}). 

PkiH  (target,  M  weapons)  =  G(cri,  <j2,  di,  d2,  h,  b2,  {fj,  j  =  1,2,...,  M})  (2) 

where  function  G  is  defined  as 

G(c 7i,a2,  d\,  d2,  6i,  b2,  {fj,  j  =  1,  2, . ,  M}) 

M 

=  E  E[F1Uu---,h)]E[F2Uu...,jk)]  (3) 

k=1 


E[F1(j1,...,jk)} 


g  A 1  (  (Pi  +  dp/k  \  * 
^1  +  ^1/  \(&i  +  +  a\) 


x  exp 


( (Eli  r«(i)A)  -  Ef=i  ^(t)2A 

2(&2  +  d2)/A; 


(Eti^(i)A)2  \ 
2((&! +  <*?)/*  +  *?) 


(4) 


/  (g+^/fc  v 

\b2  +  d\)  \(&i  +  d%)/k  +  a%) 


x  exp 


2(61  +  4)/* 


(Ef,iv(2)A-)2  N 

2((&2  +  ^2)  A  + 


(5) 


Mathematically,  F\(j\, . . . ,  jk)  is  the  product,  over  k  weapons  {j\,  ■  ■  ■ , jk},  of  all  factors 
involving  only  components  in  the  range  direction  (i.e.,  oq,  d\,  b\,  r(l))  while  F2(j  1, . . .  ,jk) 
is  the  product,  over  k  weapons  {j\, . . .  ,jk},  of  all  factors  involving  only  components  in  the 
deflection  direction  (i.e.,  cr2,  d2,  b2,  r( 2)).  Together,  equations  (2),  (3),  (4)  and  (5)  form  an 
explicit  analytical  expression  for  the  kill  probability  of  the  target,  (target,  M  weapons). 
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3  Mathematical  formulation:  multiple  weapons  against 
an  area  target  of  normally  distributed  elements 

Now  let  us  examine  the  situation  where  M  weapons  are  used  against  an  area  target 
centered  at  xtarget  —  (0,0),  consisting  of  N  discrete  elements,  normally  distributed  around 
the  center.  Let 

•  Zk  =  the  location  of  element  k  of  the  area  target,  for  k  —  1,  2, ...,  N. 

•  fj,  Xj,  and  Y  are  the  same  as  dehned  above.  They  are  respectively,  the  aiming  point, 
the  independent  error,  and  the  dependent  error  of  weapon  j. 

In  this  situation,  Zk,  the  location  of  element  k  of  the  area  target,  is  modeled  as  a  normal 
random  variable  with  zero  mean: 


We  assume  that  { Z k  —  1,2, ,  N}  are  independent  of  each  other,  and  are  independent 
of  Xj  and  Y. 

Figure  1  shows  a  sample  distribution  for  an  area  target  of  20  elements  normally  dis¬ 
tributed  with  .si  =  S2  =  4. 

To  study  the  damage  fraction  caused  by  the  M  weapons  on  the  area  target,  we  examine 
the  kill  probability  of  element  k  caused  by  the  M  weapons.  The  impact  point  of  weapon  j 
relative  to  element  k  of  the  area  target  is  given  by 

wf]  =  fj  +  Y  +  Xj  —  Zk 

=  fj  +  Y(k’eS)  +  Xj  (7) 

where  the  effective  dependent  error  of  the  M  weapons  relative  to  element  k  is  dehned  as 

y(fc,eff)  =  Y  -  Zk 
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Figure  1:  A  sample  distribution  for  an  area  target  consisting  of  20  random  elements  normally 
distributed  with  si  =  S2  =  4,  as  described  by  (6). 


Note  that  yA"’eff)  is  a  normal  random  variable  with 


zero  mean 


■y-(k,eS) 


aj  +  si 


0 


The  kill  probability  of  element  k  caused  by  the  M  weapons  (averaged  over  random  indepen¬ 
dent  errors  {Xj,  j  =  1, 2, . . . ,  M},  over  the  random  dependent  error  Y,  and  over  the  random 
element  location  Zj.)  is  given  by 


Pkin(element  k,M  weapons) 

=  G(\J of  +  s\ ,  \ja\  +  si  dx,  d2,  bu  b2,  {fj,  j  =  1,  2, . . . ,  M})  (8) 

where  function  G  is  defined  in  equations  (3),  (4)  and  (5).  Notice  that  in  the  case  of  an  area 
target  of  normally  distributed  elements,  the  kill  probability  of  element  k  has  exactly  the 
same  form  as  in  the  case  of  a  single  target  at  (0,  0)  with  the  exception  that  all  instances  of 
g\  be  replaced  by  (erf  +  sj)  and  o\  be  replaced  by  (cr|  +  s|). 


Let  Xk  be  the  Bernoulli  random  variable  indicating  whether  or  not  element  k  is  killed 
(“1”  corresponding  to  “killed”).  The  damage  fraction  (random  variable)  of  the  area  target 
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is  the  number  of  elements  killed  normalized  by  the  total  number  of  elements. 


N 


<?damage(area  target,  M  weapons)  =  —  £  Xk 

V  k=  1 


The  average  damage  fraction  has  the  expression 


<?damage(area  target,  M  weapons) 


=  E 


N 


N 

£ 

k=  1 


Xk 


N 


N 


Pkill  (element  k,M  weapons) 


k= 1 


=  G{\Jvl  +  sl  +  s2,  di,d2,  bi,  62,  {rj,  j  =  1,2,...,M})  (9) 

Expression  (9)  gives  the  exact  solution  for  the  case  of  an  area  target  of  normally  distributed 
elements. 


Figure  2:  Statistics  of  the  damage  fraction  for  the  case  of  a  single  shot  against  N  —  20 
normally  distributed  target  elements  with  no  bring  error  (ci,  cr2,  eh,  g^)  =  0  and  (61,62)  = 
(15,25).  The  radius  of  area  target  is  defined  as  Si  =  s2  (the  standard  deviation  of  the 
elements  distribution).  Left  panel:  mean  of  damage  fraction  vs  radius  of  area  target.  The 
analytical  expression  for  the  mean  is  valid  regardless  of  the  presence  or  absence  of  bring 
error.  Right  panel:  standard  deviation  of  damage  fraction  vs  radius  of  area  target.  The 
analytical  expression  for  the  standard  deviation  (equations  (11)  and  (12))  is  valid  only  in 
the  absence  of  bring  error,  which  is  true  for  the  simulations  in  this  bgure. 
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Figure  3:  Statistics  of  the  damage  fraction  for  the  case  of  a  single  shot  against  N  =  20 
normally  distributed  target  elements  with  si  =  s2  =  10,  d\  =  d2  =  0  and  (&i,62)  =  (15,25). 
The  bring  error  is  defined  as  a1  =  er2.  Left  panel:  mean  of  damage  fraction  vs  bring  error. 
The  analytical  expression  for  the  mean  is  valid  regardless  of  the  presence  or  absence  of  bring 
error.  Right  panel:  standard  deviation  of  damage  fraction  vs  bring  error.  The  analytical 
expression  for  the  standard  deviation  (equations  (11)  and  (12))  is  valid  only  in  the  absence 
of  bring  error,  which  is  false  for  the  simulations  in  this  bgure.  The  results  of  Monte  Carlo 
simulations  show  that  the  actual  standard  deviation  of  the  damage  fraction  is  larger  than 
the  value  predicted  by  applying  equations  (11)  and  (12)  with  non-zero  (oq,  cr2,  di,  d2). 

When  both  the  independent  error  (Xj,  j  =  1,2,  and  the  dependent  error  (!') 

of  the  M  weapons  are  absent,  random  variables  {xk,  k  =  1,2, .. .  ,N}  are  independent  of 
each  other.  In  this  situation,  we  can  calculate  analytically  the  standard  deviation  of  damage 
fraction  (random  variable).  The  variance  of  damage  fraction  has  the  expression 

=  var 

=  ^2  var  w  -  G)  (10) 

k= 1 

where  shorthand  notation  G  is  dehned  as 

G  =  G(\J a\  +  sf,  +  s2,  dii  ^2,  &i,  b2,  {fj,  j  =  1,  2, . . . ,  M}) 


var 


^damage (area  target,  M  weapons) 
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(11) 


Note  that  this  expression  for  the  variance  is  valid  only  in  the  absence  of  dependent  and 
independent  errors,  for  which  we  have  (ui,  cr2,  di,d2 )  =  0.  The  standard  deviation  of  damage 
fraction  is 


std 


^damage (area  target,  M  weapons) 


^-Gfl-G) 


(12) 


When  either  the  independent  errors  or  dependent  error  or  both  are  present,  the  standard 
deviation  of  damage  fraction  is  larger  than  the  value  predicted  by  applying  equations  (11) 
and  (12)  with  non-zero  (ay,  cr2l  di,  d2).  We  demonstrate  this  behavior  numerically  in  Figures 
2  and  3  below.  We  consider  the  model  problem  in  which  a  single  shot  is  bred  against  an 
area  target  of  N  =  20  normally  distributed  elements.  Monte  Carlo  simulations  are  carried 
out  with  106  runs  for  each  set  of  parameter  values  and  in  each  of  the  two  situations  below, 
yielding  accurate  numerical  results  to  compare  with  theoretical  predictions. 

In  Figure  2,  there  is  no  bring  error;  the  damage  fraction  is  affected  by  the  radius  of 
area  target  (the  standard  deviation  of  target  elements  distribution).  In  this  case,  both  the 
predicted  mean  and  the  predicted  standard  deviation  of  damage  fraction  are  valid.  As  a 
result,  the  accurate  Monte  Carlo  simulations  agree  with  the  theoretical  predictions. 

In  Figure  3,  the  radius  of  area  target  is  bxed  at  si  =  s2  =  10;  the  damage  fraction  is 
affected  by  the  bring  error  (the  total  effect  of  dependent  and  independent  errors;  for  a  single 
shot,  there  is  no  need  to  distinguish  dependent  and  independent  errors).  In  this  case,  only 
the  predicted  mean  of  damage  fraction  is  valid.  The  predicted  standard  deviation  calculated 
by  applying  equations  (11)  and  (12)  with  non-zero  (<Ji,a2,di,d2)  is  invalid  since  equations 
(11)  and  (12)  are  derived  based  on  the  assumption  of  zero  bring  error.  The  right  panel  of 
Figure  3  clearly  demonstrates  the  deviation  of  the  accurate  Monte  Carlo  simulations  from 
the  invalid  theoretical  prediction. 


4  Optimal  aiming  points  for  an  area  target 

Next  we  investigate  the  optimal  aiming  points  for  the  case  of  multiple  weapons  against 
an  area  target  of  normally  distributed  elements.  We  use  MATLAB  build-in  function  “frnin- 
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search”  [5]  together  with  formula  (9)  to  find  aimpoints  which  yields  the  largest  damage 
fraction.  In  this  study,  in  addition  to  finding  the  unconstrained  overall  optimal  aiming  posi¬ 
tions,  we  also  consider  optimizations  over  a  set  of  given  constrained  patterns  of  aiming  points. 
The  goal  is  to  find  simple  and  efficient  “empirical”  methods  for  calculating  nearly  optimal 
aiming  positions.  This  approach  greatly  simplifies  the  numerical  complexity  of  hireling  the 
optimal  aiming  points  at  the  price  of  obtaining  an  approximate  optimum.  We  consider  6 
constrained  patterns  of  aiming  positions  as  listed  below. 


Pattern  Al:  M  points  on  an  ellipse,  uniform  in  parameter  angle. 
Specifically,  the  M  aiming  points  are  mathematically  described  by 

ej  =  d1+i3^27i,  j  =  1,  2, . . . ,  M 


=  R  (^/rjcosdj,  — sin 9 


(13) 

(14) 


This  constrained  pattern  has  three  parameters:  0\ ,  R  and  r/,  over  which  we  are  going 
to  optimize  the  average  damage  fraction.  Here  6\  is  the  parameter  angle  for  weapon 
1.  R  is  the  effective  radius  of  the  ellipse,  satisfying 


R  = 


area  of  ellipse 


71 


whereas  r]  is  the  aspect  ratio  of  the  ellipse,  satisfying 

major  axis 


V  = 


minor  axis 


•  Pattern  A2:  1  point  at  center  and  (M  —  1)  points  on  an  ellipse,  uniform  in  parameter 
angle. 

One  aiming  point  is  placed  at  the  center.  The  rest  (M— 1)  aiming  points  are  distributed 
along  an  ellipse,  uniformly  in  parameter  angle  6,  as  described  by  equations  (13)  and 
(14)  where  M  is  replaced  by  (M  —  1).  This  constrained  pattern  has  three  parameters: 
6 1,  R  and  r/. 

•  Pattern  A3:  2  points  on  the  rc-axis  and  (M  —  2)  points  on  an  ellipse,  uniform  in  pa¬ 
rameter  angle. 
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Two  aiming  points  are  placed,  respectively,  at  (xa,  0)  and  (— xa,0).  The  rest  (M  —  2) 
aiming  points  are  distributed  along  an  ellipse,  uniformly  in  parameter  angle  9,  as  de¬ 
scribed  by  equations  (13)  and  (14)  where  M  is  replaced  by  (M  —  2).  This  constrained 
pattern  has  4  parameters:  9 1,  R ,  rj,  and  xa. 


Pattern  Bl:  M  points  on  an  ellipse,  uniform  in  polar  angle. 
The  M  aiming  points  are  mathematically  described  by 

0  - 1) 


=  01  + 


M 


2vr,  j  =  1,2,...,  M 


(15) 


9j  determined  by  the  condition  that  the  two  vectors 

(yrj  cos  9j,  —jn  sin  9j^J  and  (cos  0j,  sin  0j)  point  in  the  same  direction  (16) 


(Xj,Vi)  =  R  yyncosSj,  -)=sin0,j 


(17) 


This  constrained  pattern  has  three  parameters:  0i,  R  and  77. 


•  Pattern  B2:  1  point  at  center  and  (M  —  1)  points  on  an  ellipse,  uniform  in  polar  angle. 
One  aiming  point  is  placed  at  the  center.  The  rest  (M— 1)  aiming  points  are  distributed 
along  an  ellipse,  uniformly  in  polar  angle  0,  as  described  by  equations  (15),  (16)  and 
(17)  where  M  is  replaced  by  (M  —  1).  This  constrained  pattern  has  three  parameters: 
0i,  R  and  77. 


•  Pattern  B3:  2  points  on  the  rc-axis  and  (M  —  2)  points  on  an  ellipse,  uniform  in  pa¬ 
rameter  angle. 

Two  aiming  points  are  placed,  respectively,  at  (xa,0)  and  (— xa,0).  The  rest  (M  —  2) 
aiming  points  are  distributed  along  an  ellipse,  uniformly  in  polar  angle  0,  as  described 
by  equations  (15),  (16)  and  (17)  where  M  is  replaced  by  (M  —  2).  This  constrained 
pattern  has  four  parameters:  0i,  R,  77,  and  xa. 


In  numerical  simulations  below,  we  choose  the  parameter  values  as  follows. 
M  =  1  ~  12,  number  of  weapons 

(61,62)  =  (60, 100),  parameters  in  Carleton  damage  function 
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(cti,  a 2)  =  (5,5),  standard  deviation  (s.d.)  of  dependent  error  in  firing  errors 

(di,d2)  =  (5,5),  s.d.  of  independent  errors  in  firing  errors 

s  =  si  =  s2  =  15  ~  300,  radius  of  area  target  (s.d.  of  elements  distribution) 


(a) 


(b) 


(d) 


Figure  4:  Difference  in  optimal  damage  fraction  ( popt )  between  constrained  and  unconstrained 
optimizations  as  a  function  of  area  target  radius  (s)  for  various  numbers  of  weapons  (M). 
(a)  M  =  6.  (b)  M  =  7.  (c)  M  =  10.  (d)  M  =  12. 

We  compare  the  results  of  optimization  over  constrained  patterns  with  those  of  the  overall 
optimization.  Figure  4  plots  the  difference  in  optimal  (maximum)  damage  fraction  ( popt ) 
between  constrained  and  unconstrained  optimizations  as  a  function  of  area  target  radius  (s) 
for  various  numbers  of  weapons  (M).  The  difference  shows  how  accurate  it  is  to  optimize 
over  a  given  pattern.  For  M  <  6  (panel  (a)  of  Figure  4),  the  best  approximate  optimal 


13 

C-38 


damage  fraction  is  achieved  by  distribute  aiming  points  over  an  ellipse,  uniformly  in  polar 
angle  (Pattern  B1  described  above).  At  M  =  7  (panel  (b)  of  Figure  4),  the  best  approximate 
optimal  damage  fraction  is  achieved  by  placing  an  aiming  point  at  center  and  placing  the 
rest  six  aiming  points  over  an  ellipse,  uniformly  in  polar  angle  (Pattern  B2)  .  This  is  also 
true  for  M  —  8  and  M  =  9.  As  the  number  of  weapons  increases,  at  M  =  10  (panel  (c)  of 
Figure  4),  the  best  approximate  optimal  damage  fraction  is  achieved  by  placing  two  aiming 
points  on  the  x-axis  and  the  rest  eight  aiming  points  over  an  ellipse,  uniformly  in  polar  angle 
(Pattern  B3).  The  constrained  optimum  over  Pattern  B3  remains  very  accurate  at  M  =  12 
weapons  (panel  (d)  of  Figure  4). 

In  summary,  as  M  increases,  the  best  pattern  of  aiming  points  for  obtaining  approxi¬ 
mately  the  highest  damage  fraction  goes  from  Patterns  B1  to  B2  to  B3.  This  transition 
is  clearly  demonstrated  in  Figure  5  where  the  difference  in  optimal  damage  fraction  ( popt ) 
between  constrained  and  unconstrained  optimizations  is  shown  as  a  function  of  M  at  s  =  150 
(radius  of  area  target). 


Figure  5:  Difference  in  optimal  damage  fraction  ( popt )  between  constrained  and  unconstrained 
optimizations  as  a  function  of  M  at  s  =  150  (radius  of  area  traget).  In  the  horizontal 
direction,  all  points  should  have  integer  values.  To  visually  display  points  that  are  on  top  of 
each  other,  they  are  shifted  slightly  in  the  horizontal  direction  in  the  plot. 
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Figure  6:  Comparison  of  unconstrained  optimal  aiming  points  and  optimal  aiming  points 
constrained  to  Pattern  B3.  Left  panel:  two  sets  of  optimal  aiming  points  for  M  =  10, 
yielding  p{opt;  unconstrained}  =  0.772604  and  p{optj  Pattern  B3}  =  0.772398,  respectively. 
Right  panel:  two  sets  of  optimal  aiming  points  for  M  =  12,  yielding  P|op^  unconstrained}  = 
0.814997  and  pjopt  Pattern  B3}  =  0.813372,  respectively. 

Figure  6  compares  the  unconstrained  optimal  aiming  points  and  the  optimal  aiming 
points  constrained  to  Pattern  B3,  respectively  for  M  =  10  and  M  =  12  at  s  =  150.  At 
M  =  10,  the  optimal  aiming  points  of  Pattern  B3  match  almost  exactly  the  unconstrained 
aiming  points.  At  M  =  12,  the  optimal  aiming  points  of  Pattern  B3  deviate  from  the 
unconstrained  aiming  points.  Despite  the  apparent  discrepancy  between  these  two  sets 
of  optimal  aiming  points,  the  corresponding  damage  fractions  are  still  very  close  to  each 
other:  the  optimal  damage  fraction  for  Pattern  B3  is  Pjopt  Pattern  B3}  =  0-813372  while 
the  overall  optimal  damage  fraction  is  Pjopt  unconstrained}  =  0.814997.  The  difference 
between  these  two  damage  fraction  values  is  less  than  0.2%.  It  is  important  to  notice  out 
the  difference  in  computational  complexity  between  these  two  optimizations.  While  the 
constrained  optimization  over  Pattern  B3  has  4  variables,  the  unconstrained  optimization 
for  M  =  12  weapons  has  24  variables,  which  converges  much  slower  than  the  constrained 
optimization. 


15 

C-40 


The  optimal  aiming  points  constrained  to  Pattern  B1  for  M  =  6,  the  optimal  aiming 
points  constrained  to  Pattern  B2  for  M  —  7  and  the  corresponding  unconstrained  optimal 
aiming  points  are  displayed  in  Figure  7. 
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Figure  7:  Left  panel:  unconstrained  optimal  aiming  points  and  optimal  aiming  points  con¬ 
strained  to  Pattern  B1  for  M  =  6  weapons.  Right  panel:  unconstrained  optimal  aiming 
points  and  optimal  aiming  points  constrained  to  Pattern  B2  for  M  —  7  weapons. 


Figure  8:  Damage  fraction  corresponding  to  the  optimal  aiming  points.  Left  panel:  damage 
fraction  as  a  function  of  M.  Right  panel:  damage  fraction  as  a  function  of  s. 

Figure  8  plots  the  optimal  damage  fraction,  respectively,  as  a  function  of  M  for  several 
values  of  s  (left  panel),  and  as  a  function  of  s  for  several  values  of  M  (right  panel).  For  a 
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fixed  value  of  s,  the  optimal  damage  fraction  increases  with  the  number  of  weapons,  M;  for 
a  fixed  value  of  M,  the  optimal  damage  fraction  decreases  as  the  radius  (s)  of  area  target  is 
increased  (i.e. ,  damage  fraction  is  lower  for  a  larger  area  target).  Both  of  these  results  are 
reasonable  and  consistent  with  our  intuition. 

A  practical  question  regarding  resource  allocation  is  the  following.  Given  the  radius 
of  area  target  (s),  what  is  the  minimum  number  of  weapons  needed  to  achieve  a  given 
threshold  of  damage  fraction?  This  question  is  answered  in  Figure  9.  Figure  9  shows  that 
for  any  given  threshold  of  damage  fraction,  the  minimum  number  of  weapons  needed  is  an 
increasing  function  of  the  area  target  radius  (i.e.,  larger  area  target  requires  larger  number 
of  weapons),  which  again  is  reasonable  and  consistent  with  our  intuition. 


Figure  9:  Minimum  number  of  weapons  needed  for  achieving  a  given  threshold  of  damage 
fraction  vs.  radius  of  area  target.  In  the  vertical  direction,  all  flat  steps  should  have  integer 
values.  To  visually  display  steps  from  different  curves  that  are  overlap  with  each  other,  they 
are  shifted  slightly  in  the  vertical  direction  in  the  plot. 
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5  Scaling  laws  for  optimal  aiming  points  and  optimal 
damage  fraction  with  respect  to  area  target  radius 


2^1  Unconstrained  optimal  aiming  points 

2^2  Unconstrained  optimal  aiming  points 

•  Optimal  aiming  points  constrained  to  Pattern  B3 

•  Optimal  aiming  points  constrained  to  Pattern  B3 
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Figure  10:  Sets  of  optimal  aiming  points  for  M  —  10  weapons  for  several  values  of  area 
target  radius  (s).  (a)  s  =  100.  (b)  s  =  200.  (c)  s  =  300.  (d)  s  =  400. 


Finally,  we  study  how  the  optimal  aiming  points  change  with  s,  the  radius  of  area  target, 
and  explore  if  there  is  a  scaling  law  relating  sets  of  optimal  aiming  points  at  different  values 
of  s.  We  start  by  examining  the  optimal  aiming  points  for  4  different  values  of  area  target 
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radius.  The  4  panels  in  Figure  10  show  the  optimal  aiming  points  for  M  —  10  weapons, 
respectively,  for  s  =  100,  s  =  200,  s  =  300  and  s  =  400.  The  spreading  size  of  optimal  aiming 
points  increases  as  the  radius  of  area  target  (s)  is  increased.  However,  the  increase  does  not 
follow  a  simple  proportional  linear  relationship.  Figure  10  indicates  that  the  increase  in  the 
spread  size  of  optimal  aiming  points  is  less  than  linear  with  respect  to  the  area  target  radius. 
This  can  be  explained  intuitively  as  follows.  When  the  radius  of  area  target  is  increased, 
the  set  of  aiming  points  needs  to  cover  a  larger  region.  On  the  other  hand,  to  maximize  the 
damage  fraction,  the  killing  areas  associated  with  individual  weapons  also  need  to  maintain 
a  certain  degree  of  overlapping  with  each  other.  These  two  needs  contradict  each  other  and 
cannot  be  both  accommodated  simultaneously  with  a  fixed  number  of  weapons  (M)  as  the 
area  target  radius  is  increased.  Thus,  it  is  expected  that  as  the  radius  of  area  target  is 
increased,  the  spread  size  of  optimal  aiming  points  will  increase  less  than  linearly.  Here  we 
avoid  using  the  term  “radius  of  optimal  aiming  points”  because  the  distribution  of  aiming 
points  is  not  circularly  symmetric. 

For  the  purpose  of  investigating  the  spread  size  of  aiming  points  quantitatively,  we  define 
the  size  of  a  set  of  M  aiming  points  {/,  j  =  1,  2, . . . ,  M}  mathematically  as 


To  explore  how  the  size  of  optimal  aiming  points  scales  with  the  area  target  radius,  we 
plot  these  two  quantities  against  each  other  in  a  log-log  plot  in  the  left  panel  of  Figure 
11,  which  also  includes  a  fitting  function  of  the  form  Lj± p  oc  a/s.  The  log-log  plot  along 
the  fitting  function  indicate  that  the  size  of  optimal  aiming  points  (Tap)  approximately  is 
proportional  to  the  square  root  of  area  target  radius  (/s).  This  observation  suggests  that 
we  should  normalize  the  aiming  points  by  the  square  root  of  area  target  radius  and  introduce 


the  scaled  aiming  points  as 


^(scaled) 

rj 


(19) 


The  right  panel  of  Figure  11  compares  sets  of  scaled  optimal  aiming  points  of  Pattern  B3 
for  s  =  100,  s  =  200,  s  =  300  and  s  =  400.  The  comparison  demonstrates  that  not 
only  the  spread  size  of  optimal  aiming  points  scales  as  a/s,  the  distribution  of  optimal 
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Figure  11:  Comparison  of  optimal  aiming  points  for  different  values  of  area  target  radius. 
Left  panel:  spread  size  of  optimal  aiming  points  vs.  radius  of  area  target  in  a  log-log  plot. 
Right  panel:  sets  of  scaled  optimal  aiming  points  of  Pattern  B3  for  s  =  100,  s  =  200,  s  =  300 
and  s  =  400. 


aiming  points  after  scaling  is  approximately  invariant  with  respect  to  the  area  target  radius. 
Mathematically,  we  have  observed  that  approximately 


_» (scaled)  ^  .g  jnvarjan|-  respect  to  s 


(20) 


This  scaling  property  gives  us  an  even  more  efficient  way  of  calculating  optimal  aiming 
points.  We  only  need  to  calculate  the  optimal  aiming  points  for  an  area  target  of  typi¬ 
cal/representative  radius:  {fy(so),  j  =  1,  2, ... ,  M}.  We  use  so  =  150  in  our  study.  For  an 
area  target  of  radius  s,  we  simply  calculate/predict  a  set  of  nearly  optimal  aiming  points 
from  {F,(s0)}  using  the  scaling  law. 

rj{s)  =  (21) 

V  so 

We  evaluate  the  performance  of  this  efficient  method  by  examining  the  damage  fraction 
values  achieved  by  these  sets  of  nearly  optimal  aiming  points.  Specifically,  for  each  area 
target,  we  calculate  the  damage  fraction  values  corresponding  to  three  sets  of  aiming  points 


•  aiming  points  calculated  in  the  unconstrained  optimization 
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•  aiming  points  calculated  using  scaling  law  (21) 

•  all  aiming  points  =  (0,0) 


Figure  12:  Comparison  in  damage  fraction  performance  of  3  sets  of  aiming  points:  i)  aiming 
points  calculated  in  the  unconstrained  optimization,  ii)  aiming  points  calculated  using  scal¬ 
ing,  and  iii)  all  aiming  points  =  (0,0).  Left  panel:  M  =  6  weapons.  Right  panel:  M  =  12 
weapons. 

Figure  12  compares  the  damage  fraction  values  caused  by  the  3  sets  of  aiming  points 
described  above  for  M  =  6  weapons  (left  panel)  and  for  M  =  10  weapons  (right  panel).  The 
damage  fraction  achieved  by  the  set  of  nearly  optimal  aiming  points  calculated  using  scaling  is 
indistinguishable  from  that  achieved  in  the  unconstrained  optimization  (true  optimum)  while 
the  damage  fraction  corresponding  to  all  weapons  aiming  at  (0,  0)  is  much  lower.  Therefore, 
we  conclude  that  scaling  law  (21)  is  an  efficient  and  accurate  method  for  calculating  a  set  of 
nearly  optimal  aiming  points. 


6  Concluding  remarks 

We  have  studied  the  average  damage  fraction  of  an  area  target  by  multiple  weapons.  The 
area  target  was  assumed  to  consist  of  normally  distributed  elements.  Using  the  analytical 
expression  of  the  average  damage  fraction,  we  compared  various  distribution  patterns  of 
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the  aimpoints  and  gave  optimal  patterns  for  different  number  of  weapons.  Scaling  laws  for 
optimal  aimpoints  and  optimal  damage  fraction  with  respect  to  the  radius  of  the  area  target 
were  derived.  One  prospective  future  research  is  to  extend  our  current  work  to  an  area  target 
of  uniformly  distributed  elements.  Another  avenue  for  future  research  is  to  consider  an  area 
target  where  the  elements  are  assigned  different  values  and  seek  optimal  aimpoints  in  order 
to  minimize  the  total  average  surviving  value. 


7  Disclaimer 

The  authors  would  like  to  thank  TRAC-Monterey  for  supporting  this  work  and  Center 
for  Army  Analysis  (CAA)  for  bringing  this  problem  to  our  attention.  The  views  expressed 
in  this  document  are  those  of  the  authors  and  do  not  reflect  the  official  policy  or  position  of 
the  Department  of  Defense  or  the  U.S.  Government. 
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MATLAB  Code 


Damage  probability  for  multiple  weapons  against  a  unitary  target 


1  function  [ p]  =  prob_ex_sol  (M,  ra  ,  b  ,  d  ,  sig  ) 

2  Wo 

3  Wo  By  Hongyun  Wang  <hongwang@soe.ucsc.edu> 

4  Wo  University  of  California  ,  Santa  Cruz 

5  Wo  and  Hong  Zhou ,  Naval  Postgraduate  School  ,  Monterey 

e  Wo 

7  Wo  This  code  calculates  the  kill  probability  against  a  single 

8  Wo  target  at  (0,  0)  ,  by  M  weapons  with  given  aimpoints  ,  distribution 

9  Wo  of  dependent  error  ,  distribution  of  independent  errors  ,  and 
10  Wo  the  lethal  area  (Carleton  damage  function) 

n  Wo 

12  %  Input : 

13  %  M:  number  of  weapons 

14  %  ra :  M  x  2  matrix  containing  aimed  positions  of  M  weapons 

is  %  b  =  [b ( 1 )  ,  b ( 2 )  ]  :  parameters  in  Carleton  damage  function 

16  %  d  =  [d(l),  d ( 2 )  ]  :  std  of  independent  error  in  2  directions 

17  %  sig=[sig(l),  sig  (2)  ]  :  std  of  dependent  error  in  2  directions 

is  % 

19  %  Output : 

20  %  p:  kill  probability  by  M  weapons 

21  % 

22  bs=b.~2; 

23  ds=d.~2; 

24  ss=sig.  ~2; 

25  % 

26  nc=2~M— 1;  %  number  of  all  combinations 

27  cf=zeros  (nc  ,M)  ; 

28  tem=uint32  ( [  1 :  nc  ]  ’ )  ; 

29  for  k  =  1  :M, 

30  cf  ( :  ,  k)=bitget  (tem  ,  k)  ; 

31  end 

32  ka=sum(  cf  ,2)  ; 

33  % 

34  tma=(cf  *ra  ( :  ,  1 )  )  .  ~2  .  /ka  ; 

35  tmb=cf  *ra  ( :  ,  1 )  W2; 

36  Fl=sqrt  (  bs  (1) /(bs  (l)+ds  (1)  )  )Wka.*... 

37  sqrt  (  (bs (l)+ds (1) ) . /(bs (l)+ds (l)+ss (1) *ka) 

38  exp(  (tma— tmb) /(2*(bs  (l)+ds  (1)  )  )  —  tma.  /  ( 2*  (  bs  ( 1  )+ds  ( 1  )+ss  ( 1 )  *ka ) )  ); 

39  % 

40  tma=(cf  *ra  (:  ,2)  )  .  ~2  . /ka  ; 

41  tmb=cf*ra(:  ,2)  .~2; 

42  F2=sqrt  (  bs  (2) /(bs  (2)+ds  (2)  )  ).~ka.*... 

43  sqrt  (  (bs (2)+ds (2) ) . /(bs (2)+ds (2)+ss (2) *ka) 

44  exp(  (tma— tmb) /(2*(bs  (2)+ds  (2)  ) )  —  tma. /(2*  (  bs  (2)+ds  (2)+ss  (2)  *ka) )  ); 

45  % 

46  p=sum(  —  (  —  1)  Wka.  *F1.  *F2)  ; 

47  % 

48  % 
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Damage  probability  for  area  target  with  M  aimpoints  placed  along 
an  ellipse  uniformly  in  polar  angle  (EUPA) 


1 

function  [z]=p  eupa  G(M,  r, 

eta  ,  cO ,  b  ,  d ,  sig  ,  s ) 

2 

Wo 

3 

Wo  By  Hongyun  Wang  <hongwang@soe.ucsc.edu> 

4 

Wo  University  of  California  ,  Santa  Cruz 

5 

Wo  and  Hong  Zhou,  Naval  Postgraduate  School,  Monterey 

6 

Wo 

7 

Wo  This  code  calculates  the 

damage  fraction  against  an 

area 

target 

8 

Wo  of  normally  distributed 

target  elements  around  (0  , 

0) ,  by  M  weapons 

9 

Wo  with  aimpoints  placed  along  an  Ellipse  Unifromly  in 

Poar 

Angle  (EUPA)  . 

10 

Wo  The  set  of  aimpoints  is 

specified  by  3  parameters: 

major 

and  minor 

11 

Wo  axes  of  the  ellipse  ,  and 

the  off  — set  angle 

12 

Wo  Later  (in  ”p_op_eupa_G.m 

”)  ,  an  optimal  set  of  aimpoints  is  calculated 

13 

Wo  by  maximizing  the  damage 

fraction  over  the  3  parameters. 

14 

Wo 

15 

%  Input : 

16 

%  M:  number  of  weapons 

17 

%  r:  average  radius  of  ellipse  =  sqrt  (major  axis  *  minor  axis) 

18 

%  eta:  aspect  ratio  of  ellipse 

19 

%  cO  :  cO*pi  =  off-set  angle 

20 

%  b=[b(l),  b ( 2 )  ]  :  parameters  in  Carleton  damage  function 

21 

%  d=[d(l)  ,  d ( 2 )  ]  :  std  of 

independent  error  in  2  directions 

22 

%  sig=[sig(l),  sig  (2)  ]  :  std  of  dependent  error  in  2 

directions 

23 

%  s  =  [s  (1)  ,  s  (2)  ]  :  std  of  target  elements  distribution 

in  2 

directions 

24 

% 

25 

%  Output : 

26 

%  z:  damage  fraction 

27 

% 

28 

phi=cO  *  pi  +  [0  :M—  1]  ’/M*2*pi  ; 

29 

xO=cos ( phi ) / sqrt ( eta ) ; 

30 

yO=sin  (phi)*sqrt(eta); 

31 

%xO=cos ( phi ) / sqrt ( eta ) ; 

32 

%yO=sin  (phi)*sqrt(eta); 

33 

rO=sqrt  (  xO  .~2+y0 .  ~2)  ; 

34 

x=(r  *  sqrt (eta))*xO./rO ; 

35 

y=(r  /  sqrt  (eta))  *yO./rO  ; 

36 

ra  =  [x,  y  ]  ;  %  matrix  containing  aimed  positions  of  M  weapons 

37 

% 

38 

ss=sqrt  (  sig  W2+s.  ~2)  ; 

39 

[  pex]  =  prob_ex_sol  (M,  ra  ,  b  ,  cl , 

ss)  ; 

40 

z=pex  ; 

41 

% 

42 

% 
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Optimized  Damage  probability  for  area  target  with  M  aimpoints 
placed  along  an  ellipse  uniformly  in  polar  angle  (EUPA) 


1  function  [  z]=p_op_eupa_G(M,  b,  d,  sig  ,  s) 

2  Wo 

3  Wo  By  Hongyun  Wang  <hongwang@soe.ucsc.edu> 

4  Wo  University  of  California  ,  Santa  Cruz 

5  Wo  and  Hong  Zhou ,  Naval  Postgraduate  School  ,  Monterey 

e  Wo 

7  Wo  This  code  maximizes  the  damage  fraction  and  calculates  the  optimal 

8  Wo  aimpoints  when  an  area  target  of  normally  distributed  target  elements 

9  Wo  around  (0,  0)  is  attacked  by  M  weapons  with  aimpoints  placed  along 

10  Wo  an  Ellipse  Unifromly  in  Poar  Angle  (EUPA) 

n  Wo  In  the  constrained  optimization  ,  the  set  of  aimpoints  is  specified  by 

12  Wo  3  arameters  :  major  and  minor  axes  of  the  ellipse  ,  and  the  off-set  angle. 

13  Wo  The  damage  fraction  is  maximized  over  the  3  parameters. 

14  Wo 

15  %  Input : 

16  %  M:  number  of  weapons 

17  %  b=[b(l)  ,  b ( 2 )  ]  :  parameters  in  Carleton  damage  function 

is  %  d=[d(l),  d ( 2 )  ]  :  std  of  independent  error  in  2  directions 

19  %  sig=[sig(l),  sig  (2)  ]  :  std  of  dependent  error  in  2  directions 

20  %  s=[s(l)  ,s(2)  ]  :  std  of  target  elements  distribution  in  2  directions 

21  %  % 

22  %  Output : 

23  %  z:  optimal  solution  (a  structure  containing  maximum  damage  fraction  , 

24  %  optimal  aimpoints  ,  ■  ■  ■  ) 

25  % 

26  u0  =  [log  (40)  ,  log  (2)  ,  0]  ; 

27  G=@(u)  —  p_eupa_G(M,  exp(u(l)),  exp(u(2)),  u(3),  b,  d,  sig,  s); 

28  options  =  optimset  (  ’TolX  ’  ,  1  e  —  8,  ’ TolFun  ’  ,  1  e  — 16 ,  .  .  . 

29  '  Maxltcr  ’  ,2000  ,  ;  MaxFunEvals  ’  ,2000  ,  ’  Display  ’  ,  ’  off  ’ )  ; 

30  [u,  fval  ,  flag,  output]  =  fminsearch  (@(u)  G(u)  ,  uO ,  options); 

31  % 

32  pmax=—  fval  ; 

33  ap  =  [exp  (u ( l)+0 . 5  *u (2) )  ,  exp  (u  ( 1)  —  0 . 5  *u  (2) )  ,  u(3)]; 

34  rx=ap  ( 1 )  ; 

35  ry=ap  (2)  ; 

se  phi=ap  (3)  *pi  +  [0:M— 1]/M*2*  pi  ; 

37  x0=cos  (  phi )  /  rx  ; 

38  y0=sin  (  phi )  /ry  ; 

39  r0=sqrt  (  xO  .~2+y  0  W2)  ; 

40  x=rx*x0.  / r 0  ; 

41  y=ry  *  y 0 .  /  rO  ; 

42  % 

43  z=struct  ( ’M’ ,M,  ’pmax’,pmax,  ’rx’,ap(l),  ’ry’,ap(2),  ’  phi  1  ’  ,  ap  (3 )  ,  ... 

’x’  ,x,  ’  y  ,y ,  ... 

44  ’  flag  :  ,  flag  ,  ’  output  ’  ,  output )  ; 

45  % 

46  % 
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Maximum  damage  probability  of  M  weapons  against  N  targets  using 
the  optimized  EUPA  methodology 


1  Wo 

2  Wo  By  Hongyun  Wang  <hongwang@soe.ucsc.edu> 

3  Wo  University  of  California  ,  Santa  Cruz 

4  Wo  and  Hong  Zhou ,  Naval  Postgraduate  School  ,  Monterey 

5  Wo 

6  Wo  This  code  calculates  and  plots  the  maximum  damage  fraction  vs 

7  Wo  the  area  target  size. 

8  Wo  It  also  plots  the  optimal  aimpoints  for  a  given  area  target  size 

9  Wo  PROBLEM  SETUP : 

10  Wo  An  area  target  of  N  normally  distributed  target  elements  around  (0  ,0) 

n  Wo  is  attacked  by  M  weapons.  The  size  of  area  target  is  defined  as  the 

12  Wo  standard  deviation  of  the  target  element  distribution, 
is  Wo 

14  clear 

15  elf  reset 

16  % 

17  M=6;  %  Number  of  weapons 

is  N=20;  %  Number  of  target  elements 

19  % 

20 

21  b  =  [ 6 0 ,10  0];  %  parameters  in  Carleton  damage  function 

22  d=[5,5];  %  std  of  independent  error  in  2  directions 

23  sig  =  [5,5];  %  std  of  dependent  error  in  2  directions 

24  sa  =  [15:15:300]; 

25  %  A  sequence  of  values  for  the  area  target  size. 

26  %  For  each  value  ,  we  maximize  the  damage  fraction 

27  s0  =  150; 

28  %  Area  target  size  for  which  the  optimal  aimpoints  are  shown 

29  % 

30  m=si ze  ( sa  ,  2  )  ; 

31  p_eupa=zeros  ( 1  ,m)  ; 

32  for  k  =  l:m, 

33  s=sa  (k)  *  [  1  , 1  ]  ; 

34  z=p_op_eupa_G  (M,  b  ,  d  ,  s  i  g  ,  s  )  ; 

35  p_eupa  (k  )=z.pmax  ; 

36  end 

37  % 

38  figure(l) 

39  c  1  f  reset 

40  set  (  gcf  position  [50  100  560  420]) 

41  axes  (  'position  '  ,[0.15  ,0.15  ,0.75  ,0.75]) 

42  % 

43  plot(sa,  p_eupa  ,  ’  linewidth  ’  ,2 . 0  ) 

44  set  (gca  ,’ font  size  ’  ,14) 

45  axis  ([0,305,  0,  1.1]) 

46  xlabel('Area  target  size') 

47  ylabel  (  ’Maximum  damage  fraction  ’) 

48  text  (100,  1.0  ,  [num2str(M)  ,  ’  weapons  against  an  area  ... 

target  ’]  ,  ’fontsize  ’  ,14) 
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49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 


% 

figure  (2) 
elf  reset 

set  (  gcf  position  [600  200  560  420]) 
axes  (  ’position  ’  ,[0.15  ,0.15  ,0.75  ,0.80]) 

% 

s=s0  *  [1  ,  1]; 

z=p_op_eupa_G (M,  b  ,  d  ,  s i g  ,  s  )  ; 
x=z  .x  ; 

y=z  ■  y ; 

plot  (x ,  y ,  ’bh  ’  ,  ’  1  ine  width  ’  ,1.0  ,  ’markerfacecolor  ’  ,  ’y  ’  ,  ’markersize  ’  ,14) 
hold  on 

% 

ctO  =  [0:64]  / 64*2*  pi  ; 
x0=cos ( ctO  )  ; 
y0=sin  ( ctO )  ; 


plot (50*x0 ,50*y0 

5 _  1 

5  ? 

’  1  i  n  e  w  i  d  t  h  ’  ,1 

.0 

, ’color ' ,0 

.6 

*[1,  1, 

1]) 

plot (100*x0,100* 

yo, 

’  ,  ’linewidth  ’ 

,1 

.  0  ,  ’color  ’ 

,0 

-6  *  [1  , 

1,  1]) 

plot (150* xO, 150* 

yo, 

’  ,  ’  linewidth  ’ 

,1 

.0  ,  ’  color 

,0 

-6  *  [1  , 

1,  1]) 

plot (200* xO, 200* 

yo, 

’  ,  ’  linewidth  ’ 

,1 

.0  ,  ’  color 

,0 

-6  *  [1  , 

1,  1]) 

% 

set(gca, ’fontsize  ’ ,14) 
axis([  — 200, 200,  — 180,180]) 
axis  equal 

set  ( gca  ,  ’  xtick  ’  ,[—200:100:200]) 
set  ( gca  ,  ’  ytick  ’  ,[—200:100:200]) 
xlabel  (  ’Range  direction’) 
ylabel(  ’Deflection  direction’) 

text  (  —  150,  160,  ['Optimal  aimpoints  of  ’,  num2str  (M)  ,  ’ 
weapons ’] , ’fontsize  ’ ,14) 

% 

% 
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Carleton  Damage  Function  Exact  Solution  in  R 

Cardy  Moten  III 
December  27,  2016 


Purpose 

The  purpose  of  this  report  is  to  explain  the  supporting  code  for  the  equations  developed  in  Wang  et  al. 
(2016a;  2016b).  Dr.  Hongyn  Wang  and  Dr.  Hong  Zhou  originally  developed  the  code  used  in  this  report  in 
MATLAB. 


M  Weapons  firing  at  a  single  target 

Description 

The  code  in  this  section  calculates  the  damage  probability  against  a  single  target  at  (0, 0) ,  by  M  weapons 
with  given  aim  points,  distribution  of  dependent  error,  distribution  of  independent  errors,  and  the  lethal  area 
(Carleton  damage  function). 


Inputs 

•  M:  Number  of  weapons. 

•  ra:  M  x  2  matrix  containing  aimed  positions  of  the  M  weapons. 

•  =  [fo i ,  62] :  parameters  in  the  Carleton  damage  function  that  represent  the  effective  weapon  radii  in 
the  range  and  deflection  directions,  respectively. 

•  d  =  [di,  cZ2] :  standard  deviation  of  the  independent  error  in  two  directions. 

•  sig=[sig\,  sig^'.  standard  deviation  of  dependent  error  in  two  directions. 

Output 

p:  damage  probability  of  the  M  weapons  aimed  at  a  single  target  at  (0,0). 

Function  Code 

The  code  for  the  function  of  the  exact  solution  is  below.  This  code  follows  the  details,  specifically  equations 
6-8,  documented  in  Wang  et  al.  (2016a).  The  code  will  execute  the  following: 

1.  Compute  squared  values  for  parameters  6,  d ,  and  sig. 

2.  Declare  a  variable  for  the  total  number  of  hit  combinations. 

3.  Create  a  matrix  to  store  all  the  hit  combinations.  The  number  of  rows  will  be  equal  to  the  number  of 
combinations,  and  the  number  of  columns  will  equal  the  number  of  weapons  fired. 

4.  Create  a  variable  that  is  equal  in  length  to  the  number  of  hit  combinations. 

5.  For  each  row  and  column  of  the  combinations  fired  matrix: 

(a)  Compute  the  bit  value  of  the  current  index  in  the  matrix. 

(b)  For  example,  if  we  fire  two  weapons  at  a  target,  the  combinations  fired  matrix  will  contain  three 
rows  and  two  columns.  The  bit  values  for  each  row  would  be  (1,0),  (0,1),  and  (1,1)  respectively  to 
represent  a  hit  from  weapon  1,  weapon  2,  and  weapons  1  &  2  respectively. 
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6.  Compute  the  row  sum  of  the  combinations  fired  matrix. 

7.  Compute  the  impact  points  in  the  x-direction. 

8.  Compute  the  average  independent  error. 

9.  Compute  the  impact  points  in  the  y-direction. 

10.  Compute  the  average  dependent  error. 

11.  Compute  the  final  damage  probability. 

prob.ex.sol  <-  function(M,ra,b,d,sig){ 
bs  <-  b~2;  ds  <-  d~2;  ss  <-  sig~2; 
nc  <-  2~M-1  # Number  of  total  hit  combinations 

cf  <-  matrix (0,nrow=nc,ncol=M)  #Matrix  to  store  all  hit  combinations 
tmpCombinations  <-  l:nc 
for  (j  in  1:M){ 

#register  hit  and  miss  combinations 

for(k  in  1 : nrow(cf ) ){cf [k, j]  <-  as .numeric (intToBits (tmpCombinations [k] )) [j] 

> 

> 

ka  <-  rowSums(cf)  #Total  hits  per  combination  of  weapons 

#Impact  points  in  the  x-direction 
tma  <-  (cf7o*yora  [ ,  1]  )  ~2/ka 
tmb  <-  cf"/,*"/,(ra[,  1]  )~2 

#Equation  6  of  the  exact  solution  paper 

FI  =  sqrt (bs [1] / (bs [1] +ds [1] ) ) ~ka  *  sqrt ( (bs [1] +ds [1] )/ ( (bs [1] +ds [1] )+ss [1] *ka) )  * 
exp( (tma-tmb)/ (2* (bs [1] +ds [1] ) )-tma/ (2* (bs [1] +ds [1] +ss [1] *ka) ) ) 

#Impact  points  in  the  y-direction 
tma  <-  (cf7o*7«ra  [ ,  2]  )  ~2/ka 
tmb  <-  cf  (ra[,  2]  )~2 

#Equation  7  of  the  exact  solution  paper 

F2  =  sqrt (bs [2] / (bs [2] +ds [2] ) ) ~ka  *  sqrt ( (bs [2] +ds [2] )/ ( (bs [2] +ds [2] )+ss [2] *ka) )  * 
exp ( (tma-tmb) / (2* (bs [2] +ds [2] ) ) -tma/ (2* (bs [2] +ds [2] +ss [2] *ka) ) ) 

#Equation  8  of  the  exact  solution  paper 
p  <-  sum(- (-1) ~ka*Fl*F2) 
return (p) 

> 

We  will  Compute  the  single  shot  damge  probability  for  one  weapon  with  the  following  parameters  from  Wang 
et  al.  (2016a): 

.  M  =  1 

•  ra  =  [0,0] 

.  b  =  [15.4640,23.3628] 

.  d=  [5,5] 

•  sig=[30,30] 

M  <-  1;  ra  <-  matrix(c(0,0) ,ncol=2) ;  b  <-  c (15. 4640, 23. 3628) ;  d  <-  c(5,5);  sig  <-  c(30,30); 
prob . ex . sol (M , ra , b , d , sig) 
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##  [1]  0.2760992 


M  weapons  firing  at  an  area  target. 


The  code  in  this  section  calculates  the  kill  probability  where  M  weapons  are  used  against  an  area  target 
cetered  at  lit  target  =  (0,0),  consisting  of  N  discrete  elements,  normally  distributed  around  the  target  center. 
For  this  section,  we  will  implement  the  code  that  depicts  M  aiming  points  on  an  elliptical  taget,  uniform  in 
polar  angle  as  described  in  equations  15-17  of  Wang  et  al.  (2016b). 

Kill  probability  for  an  area  target 

Description 

This  code  calculates  the  damage  fraction  against  an  area  target  of  normally  distributed  target  elements 
around  (0,  0),  by  M  weapons  with  aim  points  placed  along  an  Ellipse  Uniformly  in  Polar  Angle  (EUPA).  The 
set  of  aim  points  is  specified  by  three  parameters:  weapon  effective  radius  in  the  ellipse  (?'),  aspect  ratio  of 
the  ellipse  (77) ,  and  the  off-set  angle  (cq). 


Input 

•  M :  number  of  weapons 

•  r:  average  effective  radius  of  ellipse  =  -ydnajor  axis  *  minor  axis 

•  eta(rj)\  aspect  ratio  of  ellipse 

•  Cq:  Co  *  7T  =  off-set  angle 

•  b=  [61, 62] :  parameters  in  Carleton  damage  function 

•  d  =  standard  deviation  of  independent  error  in  2  directions 

•  sig  =  [ sigi ,  sig?]:  standard  deviation  of  dependent  error  in  2  directions 

•  s  =  [si,  S2]:  standard  deviation  of  target  elements  distribution  in  2  directions 


Output 

z:  damage  fraction 


Function  Code 

To  compute  the  damage  probability  over  an  elliptical  target  area,  conduct  the  following: 

1.  Compute  a  vector  off-set  angles  for  each  weapon. 

2.  Compute  the  off-set  impact  points  scaled  by  the  aspect  ratio  in  the  x  and  y  direction. 

3.  Compute  the  radius  of  the  scaled  impact  points. 

4.  Compute  the  x  and  y  coordinates  of  the  actual  impact  points. 

5.  Generate  a  matrix  of  the  impact  points. 

6.  Compute  the  standard  deviation  of  the  weapon  error  and  target  locations. 

7.  Compute  the  average  damage  probability. 

p.eupa  <-  function(M, parameters, b,d, sig, s){ 

#Equation  15  of  the  area  target  paper 

phi  <-  parameters [3] *pi+(0 : (M-l) ) /M*2*pi  #offset  angle 

if (parameters [2] ==0){ 

parameters [2]  <-  . Machine$double . eps~ (2/3) 

> 

if  (parameters  [2]  <0)  ■[ 
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parameters [2]  <-  abs (parameters [2] ) 

> 

#Equation  16  of  the  area  target  paper 
xO  <-  cos (phi) /sqrt (parameters [2] ) 
yO  <-  sin(phi) *sqrt (parameters [2] ) 
rO  <-  sqrt (x0~2+y0~02) 

#Equation  17  of  the  area  target  paper 
x  <-  (parameters [1] *sqrt (parameters [2] )) *xO/rO 
y  <-  (parameters [1] /sqrt (parameters [2]  )) *yO/rO 

ra  =  cbind(x,y) 

ss  =  sqrt (sig~2+s~2) 

pex  <-  prob. ex. sol (M,ra,b , d, ss) 

return(-pex) 


Note,  we  are  making  the  final  probability  negative,  because  we  will  minimize  this  value  in  the  next  section. 

Optimal  aimpoints  and  Kill  probability  for  an  area  target 

The  code  below  performs  a  constrained  optimization  of  an  elliptical  target  area  that  will  produce  nearly 
optimal  aiming  points.  The  parameters  we  are  optimizing  are  the  effective  radius  (r),  aspect  ratio  (77),  and 
the  off-set  angle  (co).  According  to  Wang  et  al.,  “This  approach  greatly  simplifies  the  numerical  complexity  of 
finding  the  optimal  aiming  points  at  the  price  of  obtaining  an  approximate  optimum”  (2016b).  In  MATLAB 
Wang  et  al.  (2016a;  2016b)  used  the  fminsearch  function  for  this  optimization.  In  R,  a  similar  function  is 
found  in  the  pracma  package. 

The  code  below  executes  the  following: 

1.  Set  the  starting  conditions  for  the  optimization. 

2.  Optimize  to  find  the  effective  radius  (r),  aspect  ratio  (77),  and  off-set  angle  (co). 

3.  Store  the  maximum  probability. 

4.  Generate  a  vector  of  the  optimal  parameters  (r,  77,  Co). 

5.  Generate  the  average  radius  in  the  x  and  y  direction. 

6.  Generate  a  vector  of  off-set  angles. 

7.  Compute  the  off-set  impact  points  scaled  by  the  aspect  ratio  in  the  x  and  y  direction. 

8.  Compute  the  radius  of  the  scaled  impact  points. 

9.  Compute  the  x  and  y  coordinates  of  the  actual  impact  points. 

10.  Store  the  following  values  in  a  list: 

(a)  Number  of  weapons. 

(b)  Max  damage  probability. 

(c)  x-axis  effective  radius  (range  direction). 

(d)  y-axis  effective  radius  (deflection  direction). 

(e)  Off-set  angle. 

(f)  x-axis  aiming  points. 

(g)  y-axis  aiming  points. 

if ( ! require (pracma) ) { 

install .packages ("pracma" ) 

> 
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##  Loading  required  package :  pracma 


p.optim.eupa  <-  function(M,parameters=c(40,0,2) ,b,d,sig,s){ 
uO  <-  parameters 

opt .values  <-  fminsearch(p.eupa,c(uO[l] ,u0[2] ,u0[3] ) ,M=M,b=b,d=d,sig=sig,s=s, 

tol=le-8 ,maxiter  =  2000) 

pmax  <-  -opt . values$fval 
u  <-  opt . values$xval 

ap  <-  c (u [1] +0 . 5*u [2] ,u[l]-0.5*u[2] ,u[3] ) 
rx  <-  ap[l] 
ry  <-  ap  [2] 

phi  <-  ap [3] *pi+(0 : (M-l) ) /M*2*pi 
xO  <-  cos (phi) /rx 
yO  <-  sin(phi)/ry 
rO  <-  sqrt (x0~2+y0~2) 
x  <-  (rx*x0) /rO 
y  <-  (ry*y0)/r0 

z  <-  list(M=M,pmax=pmax,rx=ap[l] ,ry=ap[2] ,phi_l=ap[3] ,x=x,y=y) 
return(z) 


Maximum  Damage  Fraction  vs.  Target  Area  Size 

The  code  below  is  an  analysis  of  the  relationship  of  the  size  of  the  area  target  versus  the  maximum  damage 
probability.  The  initial  conditions  are  below 


•  M  =  6  weapons 

•  N  =  20  targets 
.  b  =  [60, 100] 

•  d  =  [5, 5] 

•  sig=[30,30] 

•  sa  =  sequence  from  15  to  300  in  increments  of  15  to  represent  the  size  of  the  area  target. 

The  code  and  supporting  plot  is  below.  As  shown  in  the  graphic,  as  the  size  of  the  target  increases,  the 
damage  probability  decreases  as  expected. 

M  <-  6;  N  <-20;  b  <-  c (60, 100);  d  <-  c(5,5);  sig  <-  c(5,5); 

#Generate  a  sequence  of  values  for  the  area  target  size  and 
#an  area  target  size  (sO)  to  show  the  optimal  aim  points 

sa  <-  seq(15,300,  by=  15) 

m  <-  length(sa) 
p_eupa  <-  numeric (m) 
for(k  in  l:m){ 
s  <-  sa[k]*c(l,l) 

z  <-  p.optim.eupa(M=M,b=b,d=d,sig=sig,s=s) 
p_eupa[k]  <-  z$pmax 
> 

plot(0,0,type="n" ,xlim=c(0,305) ,ylim=c(0, 1 . 1) ,xlab="Area  Target  Size" , 
ylab="Maximum  damage  fraction") 
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lines (sa,p_eupa, col="blue" , lwd=2) 

title (paste (M, "  weapons  against  an  area  target ", sep=" ") ) 


6  weapons  against  an  area  target 


Area  Target  Size 


This  final  section  of  code  plots  the  locations  of  the  optimal  aiming  points  using  the  same  initial  conditions  of 
the  previous  analysis,  however,  the  plot  will  focus  on  a  target  area  size  of  150  units2. 

sO  <-  150 
s  <-  s0*c(l,l) 

z  <-  p.optim.eupa(M=M,b=b,d=d,sig=sig,s=s) 
x  <-  z$x 
y  <-  z$y 

ctO  <-  (0:64)/64*2*pi 
xO  <-  cos(ctO) 
yO  <-  sin(ctO) 

plot(0,0,type="n" ,xlim=c(-210,210) ,ylim=c(-210,210) ,xlab="Range  direction" , 
ylab="Def lection  direction") 
lines (200*x0 , 200*y0 , type=" 1 " , lty=2) 
lines ( 150*x0 , 150*y0 , type=" 1 " , lty=2) 
lines ( 100*x0 , 100*y0 , type= " 1 " , lty=2) 
lines (50*x0 , 50*y0 , type=" 1 " , lty=2) 
points(x,y,pch=20,col="red") 

title (paste ( "Optimal  aimpoints  of  ",M,"  weapons" ,sep="")) 
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Optimal  aimpoints  of  6  weapons 


o 

o  — 

CM 


O 

O 


CM  _  ' 

1  i  i  i  i  r~ 

-200  -100  0  100  200 

Range  direction 


References 

Wang,  Hongyun,  George  Labaria,  Cardy  Moten,  and  Hong  Zhou.  2016b.  “Average  Damage  Caused  by 
Multiple  Weapons  Against  an  Area  Target  of  Normally  Distributed  Elements.”  TRAC-Monterey. 

Wang,  Hongyun,  Cardy  Moten,  Morris  Driels,  Don  Grundel,  and  Hong  Zhou.  2016a.  “Explicit  Exact  Solution 
of  Damage  Probability  for  Multiple  Weapons  Against  a  Unitary  Target.”  American  Journal  of  Operations 
Research  6  (06).  Scientific  Research  Publishing:  450. 


C-60 


Appendix  D 

Response  Surface  Methodology 


Response  surface  methodology  results 
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Precision  Guided  Munitions  Project  Notes 

Title:  Algorithm  Development  for  the  Combat  Sample  Generator  (COSAGE)  Model 
Darryl  Ahner  and  Andrew  McCarthy,  Air  Force  Institute  of  Technolgy 

Project  Goals 

The  United  States  military  has  been  using  precision-guided  artillery  rounds  (as  opposed  to 
"conventional"  rounds)  in  recent  years,  but  analytical  methods  enabling  effect  analysis  and  tradeoffs 
between  precision  and  conventional  munitions  are  lacking.  The  main  goal  of  this  project  is  to  develop  an 
algorithm  and  methodology  to  accurately  represent  the  high  explosive  (HE)  precision  munitions  using 
the  Carlton  Damage  Function. 

The  motivation  for  this  effort  is  to  accurately  model  the  desired  and  undesired  effects  of  precision 
munitions.  The  algorithm  will  be  developed  for  use  in  the  Combat  Sample  Generator  (GOSAGE)  Model  in 
use  by  the  TRADOC  Analysis  Center  -  Monterey  Naval  Postgraduate  School  and  the  Center  for  Army 
Analysis  -  Army  G-3.  The  COSAGE  model  simulates  ground  combat  between  a  large  number  of  U.S.,  Ally, 
and  Enemy  weapon  systems  over  a  48-hour  period  and  produces  engagement  results  for  U.S.  versus 
Enemy  and  or  Ally  versus  Enemy  interactions. 

The  algorithm  developed  during  this  effort  will  provide  the  Army  G-3  a  better  representation  of 
precision  munitions  effects  in  their  Joint  Integrated  Contingency  Model  (JICM)  simulation. 

Modeling  Combat  -  Initial  Setup 


Figure  1:  Range  and  Deflection  Grid 

When  speaking  about  the  grid  or  the  firing  area,  everything  is  approached  from  the  perspective  of  the 
person  firing  the  artillery  round.  To  start  with,  some  fixed  point  is  established  as  the  origin.  The  "Range" 
direction  is  determined  using  a  straight  line  from  the  weapon  to  the  target.  The  "Deflection"  direction  is 
perpendicular  to  this  line,  and  also  intersects  the  origin. 

The  entirety  of  the  analysis  is  in  metric  units  (meters).  The  standard  data  is  that  TLEs,  lethal  areas,  and 
other  measures  are  provided  with  in  metric  units. 
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For  simulation  probability  results,  a  square,  uniform  grid  is  used.  There  is  a  specified  distance  between 
each  grid  point  in  the  range  direction  and  that  same  distance  between  grid  points  in  the  deflection 
direction.  The  size  of  the  grid  is  arbitrary,  but  in  order  to  get  the  most  out  of  limited  computer  memory, 
should  be  as  small  as  possible  while  still  simulating  all  important  distribution  characteristics. 

Modeling  combat  -  Damage 


The  Carleton  Damage  Function  is  given  by  the  equation: 


P(K)  =  D0  exp  -D 


(d  d-aim ) 

2 ~b\ 


exp  —D, 


(T  raim) 

2b! 


where 

•  P(K)  is  the  probability  of  being  killed  by  that  munition  strike.  To  be  more  general,  instead  of 
saying  "kill",  "destroyed"  is  better.  This  is  because  not  all  targets  are  living.  P(K)  is  a  probability 
and  therefore  has  to  be  between  0  and  1. 

•  D0  is  the  P(K)  for  a  direct  hit.  This  is  a  function  of  both  the  target  and  the  weapon  system  and 
ranges  from  0  to  1.  For  example,  if  the  weapon  system  has  a  50%  chance  of  destroying  the 
target  with  a  direct  hit,  D0  is  equal  to  0.5 

•  d  —  daim  is  the  distance  between  the  point  in  question  and  where  the  munition  detonated  in 
the  deflection  direction 

•  r  ~  raim  is  the  same  concept,  just  in  the  range  direction 

•  b1  and  b2  represent  in  some  way  how  far  away  you  have  to  be  from  where  the  munition  lands 
to  be  safe.  To  calculate  these  two  numbers,  you  need  to  know  the  Lethal  Area  ( AL )  of  the 
munition  and  the  aspect  ratio  of  the  dangerous  area.  To  figure  out  the  aspect  ratio,  you  need  to 
know  8. 

•  8  is  the  impact  angle  of  the  munition.  This  is  the  impact  angle  relative  to  the  ground  plane;  if  the 
munition  detonates  while  pointing  straight  at  the  ground,  0  =  90°.  If  the  munition  blows  up 
parallel  to  the  ground,  0  =  0°. 


a  =  max(l  —  0.8  cos  8 , 0.3) 


bi 


la  *  Al 

x)  2n 


bz 


h. 

a 


If  the  lethal  area  (AL),  the  impact  angle  (8),  and  the  direct  hit  effectiveness  (D0),  is  known,  the  Carleton 
Damage  Function  estimates  P(K)  for  any  point  relative  to  the  explosion  point.  According  to  Klopcic 
(1990),  the  Carleton  Damage  Function  does  an  excellent  job  of  estimating  damage  done  far  away  from 
the  blast  point.  However,  that  same  paper  indicates  the  function  underestimates  damage  in  close 
proximity  as  can  be  seen  when  the  Carlton  Damage  function  is  plotted  (Figure  2). 

The  Cookie  Cutter  function  is  another  modeling  tool  used  frequently  by  the  U.S.  military.  It  assigns 
everything  in  close  proximity  to  detonation  the  same  P(K)  and  everything  farther  away  with  a  P(K)  of  0. 
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Klopcic(1990)  indicates  the  cookie  cutter  function  is  a  woefully  poor  estimator  of  outlying  damage,  but 
actually  better  than  the  Carleton  Damage  Function  for  close  proximity  damage  modeling. 
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-100 


-100 
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Figure  2:  Carlton  Damage  Function 

In  summary,  the  Carleton  Damage  Function  tends  to  underestimate  damage  in  close  proximity  to  the 
detonation  point,  but  does  an  excellent  job  approximating  damage  far  from  detonation  point,  while  the 
Cookie  Cutter  Function  is  simple  to  implement,  but  underestimates  effects  of  munitions  on  distant 
targets.  We  therefore  combine  them  taking  the  best  characteristics  from  each.  In  order  to  maximize 
how  realistic  our  simulation  is,  we  use  this  combined  function  which  is  called  the  Klopcic  Hybrid 
Function: 


Dq  if  "close" 

D0  ^exp  |d0  j  exp  ^Dq  &  ^ami)  ^  if  "not  close" 

Deciding  whether  or  not  you  are  "close"  depends  on  a  fourth  input,  the  Core  Lethal  Area  (AL0).  Using 
this  hybrid  approach  we  can  take  into  account  not  only  the  inputs,  but  also  the  how  environmental 
conditions  may  affect  those  inputs  as  indicated  in  Figure  3. 


D-4 


OUTPUTS 


CONDITIONS  INPUTS 


Terrain 

Type 

Multiplier 


P(K|Direct  Hit):  D0  — 
Core  Lethal  Area:  AL0  — 
Lethal  Area:  AL  — 

Impact  Angle:  6  — 

Deflection  distance:  d  - 
Range  distance:  r  — 


Klopcic 

Hybrid 

Damage 

Function 

1 


P k  (.Do  >  A^iAl,  9,d,r,T  LE ) 


Target  Location  Error:  TIE 


Figure  3:  Klopcic  Hybrid  Function  Inputs  and  Outputs 

The  Core  Lethal  Area  and  Lethal  Area  are  functions  of  the  weapon  system  and  the  terrain.  The  direct  hit 
effectiveness,  D0,  is  a  function  of  the  weapon  system  and  the  target  type.  This  hybrid  function  will  be  the 
basis  of  analysis  going  forward  and  can  be  used  for  several  different  estimates  of  interest.  A  three  shot 
hybrid  result  is  illustrated  in  Figure  4  where  total  probability  of  kill  as  a  function  of  location  is  depicted. 


100 


Figure  4:  Klopcic  Hybrid  Function 

As  seen  by  the  calculation,  a  closed  form  probability  of  kill  can  be  obtained  except  that  target  location 
error  exists.  This  uncertainty  changes  the  deflection  and  range  variables  from  fixed  inputs  (cf,r)  to 
random  variables  ( d+6,r+p )  where  6,p  ~  Uniform  (-TLE,  TLE)  as  seen  in  Figure  5.  To  account  for  the 
effects  of  this  target  location  error  random  variable,  several  simulation  runs  are  executed  with  the 
target  location  represented  by  uniformly  distributed  location  errors  from  0  to  the  maximum  reported 
target  location  error.  In  this  way,  the  average  effect  of  a  munition  can  be  determined.  That  effect  can 
take  the  form  of  an  object  being  in  one  of  four  states. 
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Figure  5:  Effect  of  Target  Location  Error 
Modeling  Combat  -  States  of  Combatants 

There  are  four  separate  states  that  an  object  can  exist: 

•  Unaffected  -  The  object  (whether  enemy  or  non-enemy,  mobile  or  immobile)  is  functioning  at 
the  top  of  its  game.  It  is  at  full  power. 

•  Killed  -  The  object  is  completely  destroyed  and  the  prospects  of  rebuilding  or  replenishing  the 
object  are  only  plausible  in  the  long  term. 

•  Wounded  -  An  enemy  or  non-enemy  who  is  not  killed,  but  is  slightly  damaged.  It  is  not 
functioning  at  the  highest  level  of  effectiveness,  but  is  still  armed  and  dangerous. 

•  Suppressed  -  The  object  is  not  killed,  but  is  temporarily  unable  to  defend  itself.  This  changes  as 
a  function  of  time,  but  is  considered  constant  for  the  short  time  frame  of  this  simulation.  For 
example,  if  an  enemy  is  held  down  and  cannot  return  fire,  that  enemy  is  considered  to  be 
suppressed. 

If  an  object  is  destroyed  ("killed"),  it  is  not  wounded  or  suppressed.  However,  an  enemy  can  be  both 
wounded  and  suppressed  or  either  one  alone.  Being  wounded  and  being  suppressed  are  independent  of 
one  another.  A  wounded  enemy  is  no  more  likely  to  be  suppressed  than  a  healthy  enemy.  The 
interrelation  of  these  states  is  represented  in  Figure  5. 
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Figure  5:  Interrelation  of  Targeted  Object's  States 


The  formula  for  estimating  the  probability  that  a  combatant  is  killed  is  described  earlier.  To  calculate  the 
probability  that  an  enemy  is  wounded,  the  same  logic  is  applied  as  for  this  probability  P(W),  only  the 
"Wounding"  Lethal  Area  of  the  weapon  is  larger  than  the  Lethal  Area  for  the  P(K)  calculation.  To  arrive 
at  the  "Wounding"  Lethal  Area  ,  the  lethal  area  of  the  weapon  is  multiplied  by  some  factor  dependent 
on  type  of  munition.  Because  a  wounded  object  cannot  also  be  killed,  the  Klopcic  Hybrid  analysis  for 
wounding  actually  calculates  P(W  |  Kc),  not  P(W).  This  quantity  must  be  multiplied  by  the  probability  of 
not  being  killed: 


P(W)  =  (l  -  P(/Q)  *  P(W\KC) 

The  logic  for  calculating  the  probability  of  being  suppressed  works  similarly.  The  Suppression  lethal  area 
of  the  weapon  is  the  original  lethal  area  of  the  weapon  multiplied  by  some  factor.  Generally,  this  factor 
is  larger  than  the  wounding  lethal  area  factor.  Because  suppression  is  independent  of  being  wounded, 
P(S)  is  just  multiplied  by  the  probability  of  not  being  killed.  The  equation  is: 

P(S)  =  (l-  P(K))*  P(S\KC) 


Algorithm  Development -Spreadsheet  Simulation 

KlopcicHybrid.xIsm  is  the  macro-enabled  spreadsheet  used  to  model  the  damage  done  by  a  munition  or 
group  of  munitions  while  taking  in  affect  the  uncertainty  caused  by  target  location  error.  Excel  was 
chosen  to  do  all  of  the  modeling  due  to  its  wide  availability.  It  is  limited  in  how  much  resolution  you  can 
examine  the  outcome,  but  offers  the  user  an  intuitive  set  of  inputs. 

We  describe  the  setup  of  the  simulation  in  general  here.  The  grid  setup  and  weapon  characteristics  were 
discussed  earlier  and  are  needed  to  calculate  the  Klopcic  Hybrid  function.  The  battlefield  characteristics 
are  used  to  add  a  terrain  factor  to  the  weapon  lethal  area.  Open,  level  terrain  corresponds  to  a  factor  of 
1,  while  other  terrain  (wooded,  urban,  uneven)  takes  values  less  than  1  and  shrinks  the  lethal  area  of  the 
weapon. 

The  user  sets  a  number  of  enemies  and  places  them  in  the  formation  chosen  by  the  user.  Another 
important  input  related  to  enemy  location  is  the  Target  Location  Error  (TLE).  TLE  is  not  needed  to  set  the 
grid  or  calculate  P(K)  at  each  point,  but  is  a  major  factor  in  the  analysis,  especially  when  determining 
whether  to  use  PGMs  or  conventional  rounds. 

The  value  of  TLE  quantifies  the  uncertainty  in  target  location.  In  the  real  world,  a  spotter  or  sensor  of 
some  kind  is  needed  to  determine  enemy  locations,  but  it  is  hard  if  not  impossible  to  be  exact.  For 
example,  if  the  spotter  is  fairly  certain  the  target  is  within  a  circle  of  radius  10  meters  from  where  we 
think  it  is,  the  TLE  is  10  meters.  The  target  is  equally  likely  to  be  anywhere  within  10  meters  of  the 
reported  location. 
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To  account  for  TLE  in  the  spreadsheet,  a  macro  called  PofK  goes  through  the  entire  grid  point  by  point, 
adds  up  all  the  P(K)  values  that  are  within  the  TLE  for  the  target  in  question,  and  then  divides  by  the 
number  of  grid  points  in  that  space. 

The  number  of  munitions  launched  and  how  they  are  aimed  is  the  last  input  section,  but  to  meet  the 
first  goal  of  the  project,  the  analysis  is  restrained  to  analysis  of  single  munitions.  To  describe  how 
effective  the  munition  strike  is,  several  measures  of  effectiveness  (MOEs)  are  listed.  The  expected 
number  of  kills  simply  adds  up  the  probability  of  kill  for  each  enemy  in  the  simulation.  The  expected 
number  of  kills  per  munition  divides  the  expected  number  of  kills  by  how  many  munitions  were 
launched  .The  average  P(K)  adds  up  the  P(K)  for  every  grid  point  and  then  divides  by  the  number  of  grid 
points.  Expected  number  of  wounded  and  expected  number  of  suppressed  operate  in  the  same  fashion 
as  expected  number  of  kills. 

Algorithm  Development  -  Experiment  Design 

In  order  to  measure  the  impact  of  each  variable  a  Nearly  Orthogonal  Latin  Hypercube  (NOLH)  design,  a 
space  filling  design,  is  used  to  mitigate  the  large  number  of  trials  under  various  conditions  while 
adequately  capturing  the  full  range  of  variability.  An  NOLH  ensures  the  whole  range  of  interest  of  each 
variable  is  considered,  while  having  minimal  correlation  between  input  variables.  This  approach  allows  a 
full  range  of  munition  types  to  be  considered  without  explicitly  providing  the  particular  performance 
parameters  of  munitions  as  long  as  the  parameters  are  captured  within  the  range  of  variables 
considered.  A  random  design  would  ensure  minimal  correlation,  but  wouldn't  guarantee  that  all  possible 
values  of  each  variable  are  equally  represented.  The  actual  values  used  are  given  in  Table  1. 

Table  1:  Minimum  and  Maximum  Variable  Ranges 


Input  Variable 

Minimum 

Maximum 

Deflection  Distance 

0  meters 

50  meters 

Range  Distance 

0  m 

50  m 

D0 

0.05 

1.0 

Lethal  Area  ( AL ) 

100  m2 

1500  m2 

Core  Lethal  Area 

10%  of  Al 

60%  of  Al 

Terrain  factor 

0.5 

1.0 

Impact  Angle 

30° 

CD 

O 

o 

Target  Location  Error 

2  m 

91  m 

Wound  LA  factor 

1.1 

2.5 

Suppression  LA  factor 

1.1 

2.5 

The  choice  of  minimum  and  maximum  for  each  variable  in  the  NOLH  is  very  important.  The  Core  Lethal 
Area  is  not  yet  a  measured  concept,  and  was  set  to  be  anywhere  from  10%  of  the  Lethal  Area  to  60%  of 
the  Lethal  Area.  The  terrain  coefficient  was  set  to  be  between  0.5  and  1,  and  the  wounding  and 
suppression  factors  were  set  to  be  between  1.1  and  2.5.  None  of  these  factors  had  a  meaningful  effect 
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on  the  eventual  R2  values.  Impact  Angle  can  only  vary  between  0°  and  90°,  but  any  angle  between  0°  and 
30°  provides  the  same  result  by  the  calculation.  The  assumptions  of  regression  don't  work  well  when  a 
large  subset  of  a  predictor  gives  the  same  answer,  so  the  range  chosen  was  30°  to  90°. 


That  leaves  distance  in  the  range  and  deflection  directions,  the  lethal  area,  D0  and  TLE.  The  minimum 
distance  away  is  0  meters.  You  can't  get  any  closer  than  standing  at  the  detonation  point.  The  minimum 
TLE  was  set  to  2  meters.  If  you  have  exact  knowledge  of  the  enemy's  location,  the  TLE  could  be  as  small 
as  0,  but  2  ensures  that  with  the  resolution  of  the  spreadsheet,  at  least  1  point  would  be  within  the  TLE. 

Lethal  Areas  are  classified,  so  the  estimates  are  a  relative  shot  in  the  dark.  In  the  example  that  Major 
Moten  gave  us,  the  lethal  area  was  2,270  ft2  (211  m2).  Based  on  this,  the  minimum  value  was  set  at  100 
m2. 

In  terms  of  the  maximum  value  for  TLE  and  lethal  area,  we  discussed  with  Mr.  Robert  Lillard  of  the  Fires 
Center  of  Excellence  whether  our  choices  of  maximum  TLE  and  maximum  lethal  areas  were 
representative  of  real  values.  A  maximum  lethal  area  of  1500m2  and  a  maximum  TLE  of  25m  (for  the 
most  precise  weapon  system,  Excalibur)  were  given  a  reasonable  stamp  of  approval.  It  was  then  at  our 
discretion  to  determine  how  much  distance  we  want  to  analyze  and  what  values  of  D0  to  use. 


Variable 

Min 

Max 

Deflection  (d) 

Om 

50  m 

Range  (r) 

Om 

50  m 

Do 

0.05 

1 

CLA 

0.1 

0.6 

LA 

100  m2 

1000  m2 

Tf  actor 

0.5 

1 

ImpAng 

30° 

90° 

TLE 

2m 

91  m 

Table  2:  Initial  Experimental  Range 


There  are  three  types  of  tube-launched  artillery  rounds  in  question:  Excalibur,  PTK,  and  conventional 
rounds.  Deciding  which  one  to  fire  is  mostly  a  function  of  the  TLE.  If  the  TLE  is  25  meters  or  less,  two 
Excalibur  rounds  are  fired  at  the  target.  These  currently  cost  around  $5,000  per  round,  but  presumably 
may  become  more  expensive  in  the  near  future  ($15,000-$20,000  or  so).  If  the  TLE  is  greater  than  25 
meters  but  less  than  100  meters,  PTK  rounds  are  fired.  These  are  not  as  precise  as  Excalibur,  but  are  still 
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considered  close  to  "precision-guided  munitions".  How  many  are  fired  depends  on  the  target,  but  these 
are  generally  fired  in  small  groups. 

Any  TLE  above  100  meters  will  be  fired  upon  using  conventional  rounds,  generally  in  groups  of  30-40 
rounds  per  attack.  "Dumb"  rounds  cost  around  $1000  per  round.  After  this  initial  firing,  there  is  a  re¬ 
assessment  period  and  the  decision  is  made  whether  to  attack  again.  Whatever  kind  of  round  is  being 
fired,  the  target  is  called  a  unit.  A  unit  could  be  anything  from  a  small,  single  stationary  object  to  several 
vehicles  50  to  75  meters  apart. 

The  simulation  of  these  weapon  systems  can  be  three-dimensional,  but  the  algorithm  should  not  be.  The 
Carleton  Damage  Function  is  not  designed  to  be  a  three-dimensional  equation.  Also,  the  Excalibur 
weapon  is  not  designed  to  perform  area  fires,  where  many  rounds  are  fired  at  a  group  of  targets  in  a 
large  area.  Conventional  rounds  are  best  suited  for  this  purpose. 

PK  =fi(d,  r,  D0,Al,  Alo , T,  0,  TLE) 

Equation  1 

Initial  experimentation  is  conducted  to  determine  the  overall  fit  of  a  regression  model  to  a  full  range  of 
factors.  If  the  model  is  sufficient  at  this  point  no  further  model  would  be  required. 


a  Actual  by  Predicted  Plot 


A  Summary  of  Fit 

RSquare  0.720011 

RSquare  Adj  0.6619 

Root  Mean  Square  Error  0.024127 

Mean  of  Response  0.031862 

Observations  (or  Sum  Wgts)  257 


Analysis  of  Variance 

Source 

DF 

Sum  of 
Squares 

Mean  Square 

F  Ratio 

Model 

44 

0.31734018 

0.007212 

12.3903 

Error 

212 

0.12340328 

0.000582 

Prob  >  F 

C.  Total 

256 

0.44074345 

<.0001' 

Figure  6:  Initial  Model 
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Algorithm  Development  -  Experiment  Results 

A  NOLH  was  constructed  using  the  ranges  in  Table  2  for  the  initial  experiment.  The  initial  model  fitted  a 
response  service  that  allowed  polynomial  of  degree  terms  and  two-way  interaction  terms  to  be 
considered  by  a  step-wise  regression  and  is  a  function  of  terms  shown  in  Equation  1.  This  resulted  in  a 
marginal  fit  with  an  adjusted  R2  of  .66  which  is  less  than  expected  and  is  shown  in  Figure  6.  Additionally, 
the  data  displayed  a  nonlinear  relationship  with  the  fitted  function  indicating  a  transformation  of  the 
data  may  be  required.  A  square  root  transformation  was  applied  to  the  dependent  variable  P(k)  to  seek 
to  correct  this  nonlinear  relationship. 


J  Actual  by  Predicted  Plot 


sqrtPK  Predicted  P<.0001  RSq=0.83  RMSE=0.0442 


-J  Summary  of  Fit 

RSquare 

0.829603 

RSquare  Adj 

0.794238 

Root  Mean  Square  Error 

0.044186 

Mean  of  Response 

0.149701 

Observations  (or  Sum  Wgts) 

257 

J  Analysis  of  Variance 

Sum  of 


Source  DF  Squares  Mean  Square  F  Ratio 

Model  44  2.0152060  0.045800  23.4581 

Error  212  0.4139140  0.001952  Prob  >  F 

C.  Total  256  2.4291199  <.0001* 


Figure  7:  Initial  Transformed  Model  SQRT(P(k)) 


The  transformed  P(k)  regression  is  shown  in  Figure  7.  Despite  the  increased  adjusted  R2  value  of  .79, 
there  are  concerns  over  whether  this  was  a  good  fit  across  all  variable  values  or  whether  it  was  just  a 
good  fit  for  very  low  P(K)  values  far  from  the  target  and  mediocre  for  close  proximity  points  which  have 
high  P(K)  as  shown  for  the  circled  values  in  Figure  7.  To  mitigate  this,  models  were  developed  for 
different  maximum  values  of  TLE,  Lethal  Area,  and  distance  from  the  origin  (d  and  r).  Specifically 
different  models  were  developed  for  maximum  TLEs  of  15m,  50m,  and  91m,  maximum  Lethal  Areas  of 
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500m2,  1000m2,  and  1500m2,  and  maximum  d  and  r  values  of  10m,  15m,  20m,  25m,  and  50m.  In  all, 
models  were  developed  for  45  different  testing  configurations.  The  goal  was  to  examine  how  R2  is 
affected  by  simply  changing  the  range  of  several  key  input  variables  using  R2  as  the  dependent  variable. 
The  result  of  this  analysis  was  to  verify  that  models  with  large  TLEs,  small  Lethal  Areas,  and  large  values 
of  d  and  r  produce  poor  predictive  models.  While  not  surprising,  it  highlights  that  in  order  to  develop  a 
realistic  model,  the  ranges  of  input  variables  had  to  be  realistic  and  these  parameter  settings  from  the 
initial  model  do  not  add  resolution  to  the  model. 

Using  this  insight,  separate  models  were  developed  for  Excaliber,  PGM,  and  conventional  munitions 
using  values  less  than  25  meters,  100  meters,  and  1000  meters  respectively.  Since  from  the  original 
transformed  model  and  further  experimentation  it  was  known  that  model  fit  may  be  effected  by 
nonlinearity  in  some  key  variables,  classification  and  regression  trees  were  used  to  determine  which 
variable  created  the  most  variability  in  the  model.  As  shown  for  the  dataset  with  TLE  <  25,  there  is  a 
natural  break  at  the  deflection  of  12.6  meters.  This  process  was  subsequently  used  to  divide  the  data  as 
shown  in  Table  3. 


a  »  Partition  for  sqrtPK 
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;  ,  *•  * 
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»  •  ,  •  ,  • 
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.  •  .  •  .  *  * 

Def>=12.6 

Def<12.6 

All  Rows 


Split  j  |  Prune 


Number 

RSquare 

RMSE 

N  of  Splits 

AICc 

0.362 

0.1456852 

257  1 

-254.69 

"  All  Rows 

Count  257  LogWorth  Difference 

Mean  0,3742079  39.467181  0.21935 

Std  Dev  0.1827073 


’  Def> 

Count 

Mean 

Std  Dm 

=12j6 

128 

0.2641086 

0.140515 

l’Def<12j6 

Count  129 

Mean  0.4834536 
Std  Dev  0.1517365 

>  Candidates 

>  Candidates 

Figure  8:  Classification  and  Regression  T ree  for  TLE<25 
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This  process  is  followed  to  develop  models  for  both  Excaliber  and  PGK  munitions  as  shown  in  Figure  9 
(See  Appendices  for  equations)  .  This  methodology  creates  the  4  models  shown  in  Figure  9  with  model 
goodness-of-fits  that  are  both  traceable  and  defendable  for  munition  effects.  The  Fit  adjusted  R2  values 
are  shown  with  the  highest  R2  values  for  the  models  with  rounds  landing  closest  to  the  target  regardless 
of  the  lethal  area  which  is  the  area  of  most  importance.  The  Fit  adjusted  R2  values  are  for  the  modeling 
of  the  square  root  of  P(k).  Validation  R2  values  are  for  simulations  of  actual  P(k)  after  the  square  root  is 
transformed  back  to  P(k).  These  validation  runs  used  500  randomly  selected  factor  vectors.  Again,  the 
model  performs  extremely  well  for  the  area  of  most  concern  of  higher  P(k)  values  close  to  the  target. 


Model 

Conditions 

Model  1A 

Dand  R<12.6,  TLE  <  25 

Model  IB 

D  or  R  >  12.6,  TLE  <  25 

Model  2A 

Dand  R  <  25,  100  >  TLE  >  25 

Model  2B 

D  or  R  >  25,  100  >  TLE  >  25 

Table  3:  Data  Breakpoints  for  Models 


Fit  (sqrt(Pk))  R2=.95  R2=.89  R2=.99  R2-82 

Validation  (Pk)  R2=.95  R2=.77  R2=.96  R2=.67 

Figure  9:  Final  Model  Goodness-of-Fit  and  Validation  Results 

Adjustment  for  Risk  Averse  Modeling 

While  the  preceding  models  capture  well  the  average  P(k)  effect  of  a  single  mention  round,  planners 
sometimes  want  to  be  risk  averse  do  to  civilians  or  proximity  of  sensitive  objects  to  impact.  As  with  any 
good  regression  model,  the  preceding  models  had  errors  in  prediction  that  were  verified  to  be  normally 
distributed  .  For  risk  averse  modeling,  we  use  this  normal  distribution  of  errors  to  adjust  our  model. 
Based  on  the  Model  in  use,  a  constant  based  on  error  distribution  can  be  added  to  the  model  before  it  is 
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squared  and  then  both  the  model  function  and  added  constant  to  predict  a  risk  based  P(k)  as  shown  in 
Equation  2. 


=  f2(d,r,D0,AL,  Aw  ,T,0,  TLE)  +  C 
PK  =  [f2  (d,  r,  D0 ,  Al ,  Aw , T,  e: TLE)  +  C]2 

Equation  2 

The  constant  C  is  determined  through  a  percentile  plot  of  the  error  terms  shown  in  Figure  10. 
Dependent  on  the  level  of  risk  aversion,  C  values  for  each  model  can  be  obtained.  For  Example,  values 
for  C  are  given  for  the  75th,  90th,  and  99th  percentile  in  Table  4. 


Model  1 A  Model  IB  Model  2A  Model  2B 

4 1  Residual  sqrtPK  £  -  Residual  sqrtPK  2  A  *  Residual  sqrtPK  2  d|g  Residual  sqrtPK 


Figure  10:  Error  for  Each  Model 


Model 

75%  Adj 

90%  Adj 

99%  Adj 

1A 

0.014 

0.032 

0.10 

IB 

0.04 

0.077 

0.15 

2A 

0.005 

0.011 

0.04 

2B 

0.022 

0.041 

0.10 

Table  4:  Error  for  selected  percentiles  for  each  model 


Application  of  Models 

To  apply  the  models ,  traceable  inputs  for  each  scenario  must  be  known.  Given  this  information,  each 
round  is  assumed  to  act  independently.  This  model  is  typical  in  this  type  of  modeling  and  does  not 
account  for  a  situation  where  one  round  may  wound  an  individual  and  that  same  individual  is 
"wounded".  The  number  of  rounds  fired  in  each  situation  is  left  to  the  modeler;  however,  this  modeling 
approach  allows  the  modeler  to  begin  to  approach  the  question  of  how  many  rounds  are  enough  to 
achieve  the  desired  effect. 
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Conclusion 


This  effort  developed  an  algorithm  and  methodology  to  accurately  represent  the  high  explosive  (HE) 
precision  munitions  using  the  Carlton  Damage  Function  for  use  with  the  COSAGE  model  to  accurately 
model  the  desired  and  undesired  effects  of  precision  munitions.  This  effort  should  enable  the  COSAGE 
model  to  better  simulate  precision  munition  weapon  effects  in  ground  combat  between  a  large  number 
of  U.S.,  Ally,  and  Enemy  weapon  systems  over  a  48-hour  period  and  produce  engagement  results  . 

Similar  models  can  be  developed  for  probability  of  wounding  and  probability  of  suppression. 
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APPENDIX  1  -  Excaliber  Models 


Trial  55 

Variable 

Min 

Max 

Deflection 

0  m 

25  m 

Range 

0  m 

25  m 

Do 

0.25 

1 

Lethal  Area 

100  m2 

1500  m2 

Core  Lethal  Area 

0.1 

0.6 

Terrain  Factor 

0.5 

1 

Impact  Angle 

30 

90 

TLE 

2 

25 

APPENDIX  1A-  Excaliber  Model  1A 
SQRT(P(k)  = 

0.3 15 1885-0.012095  *DEF- 

0. 006095  *RAN+0.2659346*D0+0. 0215783  *CLA+0.0002081*LA+0.1779652*TFACTOR+0. 0008161*  I  F(  I M 

PANG<30,30,IMPANG)-0.007863*TLE-0.000866*(DEF-5.73838)A2-0.000329*(RAN-5.77882)A2- 

0.014571  *(DEF-5. 73838)  *(D0-0.62775)-0. 009837  *(RAN-5.77882)*(D0-0.62775)-0.214067*(D0- 

0.62775)  A2+0. 0000047582*  (DEF-5. 73838)  *(LA-804.634)+0. 0000035 178*  (RAN-5. 77882)*  (LA- 

804. 634)+0. 0002606*  (DO-0.62775)  *(LA-804.634)-0.0000001878*(  LA-804. 634)  A2+0.2768908*(  DO- 

0.62775)*(TFACTOR-0.75083)-0.0000738*(LA-804.634)*(TFACTOR-0.75083)-0.220381*(TFACTOR- 

0.75083)A2+0.0001185*(DEF-5.73838)*(IF(IMPANG<30,30,IMPANG)-60.0127)-0.00008652*(RAN- 

5.77882)*(IF(IMPANG<30,30,IMPANG)-60.0127)+0.0014057*(D0- 

0.62775)*(IF(IMPANG<30,30,IMPANG)-60.0127)-0.000022*(IF(IMPANG<30,30,IMPANG)- 

60.0127)  A2+0.0012364*(DEF-5.73838)*  (TLE-13. 535)+0. 0004483  *(RAN-5. 77882)*  (TLE-13.535)- 

0. 008522  *  (DO-0. 62775)  *(TLE-13. 535)-0. 000226*  (TLE-13. 535) 

For  30  degrees  <=  alpha  <=  90 

APPENDIX  IB-  Excaliber  Model  IB 

SQRT(P(k)  = 

0. 1821638-0. 007924*  DEF- 

0.00459*RAN+0.000169*LA+0.160827*TFACTOR+0.0003699*IF(IMPANG<30,30,IMPANG)+0.0023676*T 
LE+0.000081642*(DEF-26.9887)A2+0.000089675*(DEF-26.9887)*(RAN-26.9609)+0.000024447*(RAN- 
26.9609)  A2-0. 00128*  (RAN-26. 9609)  *(  D0-0.623 14)-0. 000003489*  (DEF-26.9887)*(  LA-796. 842)- 
0. 000002  *(RAN-26. 9609)*  (LA-796.842)+0.00006377*(  D0-0.623 14)*  (LA-796.842)-0.00000006079*(  LA- 
796.842)  A2-0. 003067  *(DEF-26. 9887)*  (TFACTOR-0. 7495  )+0.0001227*(  LA-796.842)  *(TFACTOR-0. 7495)- 
0.00003078*(RAN-26.9609)*(IF(  IM  PANG<30, 30,  IM  PANG)-60. 0044)-0. 00004532*(  RAN-26. 9609)*(TLE- 
13. 4792)+0. 0055151*  (D0-0.623 14)*  (TLE-13. 4792)-0.000003057*(  LA-796.842)*  (TLE-13. 4792) 
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APPENDIX  2 -PGK  Models 


Trial  57 

Variable 

Min 

Max 

Deflection 

0  m 

100  m 

Range 

0  m 

100  m 

Do 

0.25 

1 

Lethal  Area 

100  m2 

1500  m2 

Core  Lethal  Area 

0.1 

0.6 

Terrain  Factor 

0.5 

1 

Impact  Angle 

30 

90 

TLE 

2 

25 

APPENDIX  2 A-  PGK  Model  2A 
SQRT(P(k)  = 

0.2 118496-0.0005  6*  DEF- 

0. 000298*  RAN+0. 0247488*  D0+0. 0001617*  LA+0.1465859*TFACTOR+0. 0001839*  IF(IMPANG<30,30,IMP 
ANG)-0.003926*TLE-0.00003766*(DEF-9.1624)A2-0.00002676*(RAN-8.86918)A2-0.000000894*(DEF- 

9.1624) *(LA-790.25)-0.0000005558*(RAN-8.86918)*(LA-790.25)+0.0000192*(D0-0.61069)*(LA-790.25)- 
0.00000007661*(LA-790.25)A2+0.0000823*(LA-790.25)*(TFACTOR-0.74579)-0.103632*(TFACTOR- 
0.74579)A2+0.00000024725*(LA-790.25)*(IF(IMPANG<30,30,IMPANG)-59.9399)+0.0000414*(DEF- 

9.1624) *(TLE-58.6302)+0.000044243*(RAN-8.86918)*(TLE-58.6302)-0.00068*(D0-0.61069)*(TLE- 
58. 6302)-0. 000002 11*  (LA-790. 25)  *(TLE-58.6302)-0. 001881  *(TFACTOR-0.74579)*(TLE-58.6302)- 
0.00001285*(IF(IMPANG<30,30,IMPANG)-59.9399)*(TLE-58.6302)+0.000058223*(TLE-58.6302)A2 


APPENDIX  2B-  PGK  Model  2B 
SQRT(P(k)  = 

0.057898 1-0.001659*  DEF- 

0.001348*RAN+0.0000672*LA+0.0735056*TFACTOR+0.0011076*TLE+0.0000294*(DEF-53.0857)*(RAN- 
52.9462)-0. 000001 108*  (DEF-53.0857)*(  LA-805. 155)-0. 0000007203  *(RAN-52.9462)*(  LA-805. 155) 
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Appendix  F 
Glossary 


ARO 

CAPTTIM 

IGT 

TRAC 

WCST 


Army  Research  Office 

Cognitive  Alighnment  With  Performance  Targeted 
Training  Intervention  Model 
Iowa  Gambling  Task 

Training  and  Doctrine  Command  Analysis  Center 
Wisconsin  Card  Sorting  Test 
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